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EXCESS PROPERTIES OF H, -D, LIQUID MIXTURES* 


M. Lambert 
Faculty of Sciences, University of Brussels, Brussels, Belgium 
(Received April 18, 1960) 


The heat of mixing (h,) and excess volume (v,) 
of H, - D, liquid mixtures have been measured 
at 20.4°K. The method used is similar to that 
described in detail by Jeener,' and a complete 
description of the apparatus will soon be sub- 
mitted for publication. 

The gaseous H, and D, are condensed in a 
calorimeter on each side of a tight diaphragm 
which may be broken from outside. The mixing 
takes place at constant pressure (1.5 atm) and 
adiabatically. The calorimeter temperature is 
measured by means of a platinum resistance 
thermometer (500 ohms at room temperature, 

12 ohms at about 20°K) and a Mueller bridge. 

A temperature change of one degree gives rise 
toa resistance variation of nearly 0.5 ohm, and 
changes of 0.0002 degree may so be observed. 
The thermometer is calibrated with reference 
tothe vapor pressure of H,. A heating coil per- 
mits the heat capacity of the calorimeter to be 
determined. The heat generated by breaking the 
diaphragm and by stirring is small compared to 
the heat of mixing; corrections for these effects 
are made. 

The hydrogen was obtained from the Society 
“L’Air liquide,” Liége, Belgium, and was re- 
ported to be 99.99 % pure. 

The deuterium was obtained by electrolysis of 
heavy water. Mass spectrometer analysis showed 
it to contain about 1% HD. 

The measurements were made with liquids in 
ortho-para equilibrium at 20°K (99.8% p-H,, 

%% o-D,). Conversion was achieved in the 
liquid phase according to a method described 


elsewhere.* The composition was determined in 
the gas phase by thermal conductivity analysis. 
The results are listed in Table I. 
Newman and Jackson* obtained for the excess 
free energy of this system the value 


a 1.6+0.1 cal/mole at 22°K. 


Although these measurements are not at the same 
temperature as ours, we may safely conclude that 
the excess entropy s, > 0. 

The theory of Prigogine and co-workers for 
isotopic mixtures*’® predicts for this system, at 
20.3°K and x =0.5: 


&,= 1.85 cal/mole, Ts, =2.2 cal/mole, 
h, =4.05 cal/mole, v, = 0.21 cm*/mole. 
While the excess enthalpy is in satisfactory 


agreement, the excess volume has the opposite 


Table I. Summary of results. 





x(H,) v,(em?/mole) he (cal/mole)* hp/x(H,)x(D,) 





0.300 -0.14 2.56 12.2 
0.436 -0.24 2.93 11.9 
0.229 2.03 11.5 
0.559 2.92 11.8 
0.688 -0.16 2.56 11.9 
0.804 1.79 11.4 





“The excess enthalpy may be represented by the 
following equation: he= ax(H,)x(D,), with a =11.8 
cal/mole. The discrepancies do not exceed 5%. 
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sign. In the theoretical predictions the sign of 
the excess volume changes according to the mag- 
nitude of the quantum effects: a contraction is 
expected for large effects (He* - He* system) 
and a dilatation for small effects. It is rather 
difficult to ensure which is exactly the position 
of the H, -D, system in this problem. It is also 
necessary to bear in mind that the theory was 
developed for monatomic molecules, which is 
not the case here. 

Kerr® observed for the similar D, -T, system 
a positive excess volume. 

Measurements of hz and vg were also made on 
the o-H,—p-H, system. For equimolar mixtures 
and at 20.4°K, we found 


h, =0.4 cal/mole; 
a slight dilatation of the order of 0.1% was ob- 
served. The value of h, is of the same order of 


magnitude as the value calculated by Friedman’ 
from vapor pressure data®: 


h, =0.7 cal/mole. 


More extensive measurements over various 


temperatures are in progress, which will give 
information on the excess specific heat. Also 
measurements on Ne”° - Ne” mixtures will be 

made for the same properties. 

We wish to thank Professor Prigogine for con- 
stant interest on this work. We are also in- 
debted to J. Jeener, V. Mathot, and A. Bellemans 
for many interesting discussions. 





"This work has been supported by grants from the 
“Institut pour 1’ Encouragement de la Recherche 
Scientifique dans 1’Industrie et l’Agriculture” and 
“Le Centre National Belge de Chimie Physique 
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EXPERIMENTAL TWO-BEAM EXCITATION OF ELECTRON OSCILLATIONS 
IN A PLASMA WITHOUT SHEATHS 


M. J. Kofoid 
Boeing Scientific Research Laboratories, Seattle, Washington 
(Received April 15, 1960) 


Coherent standing waves of longitudinal elec- 
tron oscillations have been excited in a plasma 
by two independent oppositely-directed electron 
beams, with no sheaths on the electrodes from 
which the electron streams entered the plasma. 
The axes of the two beams coincided. 

The absence of sheaths and the independence of 
the beams made this a very different experiment 
than that of Looney and Brown.' 

In the initial tests, identical beams of 65-ev 
electrons entered the plasma through 1.0-mm 
diameter apertures in plane electrodes. These 


electrodes were at the potential of the plasma, or 


slightly positive, and therefore were without 
sheaths. Both beam electrodes were at the same 
potential. The apertures faced each other across 
2.0 cm of plasma whose electron density corre- 
sponded to a Langmuir frequency of about 550 

Mc/sec. The electron density of each beam was 
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of the order of 10% of the density of this back- 
ground plasma. The tests were made in argon at 
11 » pressure. Detection was by means of a 
0.20-mm diameter movable probe located 2 mm 
out from the visible edge of the interpenetrating 
beams. Strong, reasonably sinusoidal standing 
waves, with two loops between beam electrodes, 
were found at 550 Mc/sec when (and only when) 
both beams were present. The standing waves 
could not be produced without the background 
plasma between the beam electrodes. With the 
distance between electrodes taken as a wave- 
length, the phase velocity was 1.10 10° cm/sec; 
the electron velocity of either beam was 0.48x10 
cm/sec. 

With the two-loop mode of operation, the oscil- 
lation frequency spectrum was as shown in Fig. |. 

A single-loop standing wave was established 

















with a beam energy of 74 ev. With either one or 
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FIG. 1. Oscillation frequency spectrum as por- 


trayed on the screen of a spectrum analyzer. Center 
frequency is 550 Mc/sec. 


two loops, the standing-wave amplitude fell grad- 
lly to zero when the voltage was increased or 
decreased about 3 volts from the stated values. 
The principal elements of the experimental 
arrangement are shown in Fig. 2. In each of the 
identical guns the source of electrons was a 200- 
ma discharge between the heated, coated cathode 
Cand the anode A,; electrons extracted through 
a1.0-mm diameter aperture from this discharge 
were accelerated and focused by the electric 
field between electrodes A, and A,. The gun arc 
discharges were operated at 25 » pressure. The 
current in each beam as it left electrode A, was 
aout 0.5 ma. The background plasma between 
the beam electrodes was produced by running an 
arc discharge of about 300 ma total between two 
identical wire ring anodes A’ (which were con- 
nected directly to each other) and two heated, 





These cathodes were located diametrically oppo- 
site to each other. At each plasma discharge 
tathode C’, at a distance of a few mm from it, a 
shield P was located so as to prevent electrons 
from traveling without collision to the region 
penetrated by the beams. 

The possibility was realized that electrons from 
one gun could enter the aperture of the second 
gun to be reversed to full velocity by the accele- 
tating field of the second. Thus, a reverse beam 
might be created for each initial beam. It was 








coated cathodes C’, only one of which is indicated. 
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FIG. 2. Test apparatus arrangements. The A, 
elements have been called the “beam electrodes.” 
Pumping line was throttled to control pressure in main 
chamber. 


reasoned that perhaps “two-beam operation” was 
being established independently by each initial 
beam. However, when the accelerating voltage 
was left on both guns but the arc discharge was 
extinguished in one gun, it was found that no elec- 
tron oscillations could be excited by the operating 
beam as it penetrated the plasma. It was then 
concluded that not only the electrons reflected by 
the field, but also the fewer electrons reflected 
by impact on electrodes, were playing no im- 
portant role in the experiment. 

The single-beam experiments of Looney and 
Brown’ with equal sheaths and of Gordon’ with 
unequal sheaths on the gun and reflector elec- 
trodes were verified with the apparatus used for 
the two-beam experiment with no sheaths. The 
basic elements were constructed essentially iden- 
tical to those of Looney and Brown in order that 
their work could be checked directly. 

Credit belongs to Dr. J. E. Drummond for point- 
ing out the need of experiments of the general na- 
ture of that performed and for numerous discus- 
sions of the phenomena and theory involved. 
Valued discussions with Professor K. G. Emeleus 
and Professor R. Geballe are gratefully acknowl- 
edged. Dr. S. Sobottka aided through helpful con- 
versations and apparatus improvements. C. W. 
Shirley and H. V. Cleva gave continual assistance 
in the experimental work. 





'D. H. Looney and S. C. Brown, Phys. Rev. 93, 
965 (1954). 

?E. I. Gordon, Ph.D. thesis, Massachusetts Insti- 
tute of Technology, 1957 (unpublished). 
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MEASUREMENT OF THE POLARIZATION OF SYNCHROTRON RADIATION* 


Peter Joost 
Cornell University, Ithaca, New York 
(Received May 9, 1960) 


The radiation emitted by electrons accelerated 
in a synchrotron has been investigated theoreti- 
cally’»? and its spectrum and angular distribution 
have been measured.* It has an important astro- 
physical aspect since there is strong evidence 
that the light of the Crab Nebula has its origin in 
high-energy electrons moving through a strong 
magnetic field. Measurements on the continuous 
spectrum of the Crab Nebula have shown the light 
to be partly linearly polarized.* However, at this 
time one assumed that synchrotron radiation was 
linearly polarized. Recently Westfold® has shown 
for electrons moving at an angle a to a magnetic 
field and Bethe® for electrons in a synchrotron 
that the radiation is elliptically polarized, the 
ratio between major and minor axes depending 
on the angle ~ between direction of observation 
and the synchrotron plane. This makes it desir- 
able to measure the polarization of synchrotron 
radiation from a terrestrial source. 

The amplitude of light whose plane of polariza- 
tion lies in the synchrotron plane is® 


= = fw os 
Je° "tg * 2/3 (5, , 


and of that with plane of polarization perpendic- 
ular to it, 


— ‘a 
J ai *¥®/3(555*) 


Here x =(¢7+y~*)” with y =1/(1 - v?/c?)”, where 
v is the velocity of the electrons, c that of light, 
the K’s are modified Bessel functions of the sec- 
ond kind, w is the angular frequency of the ob- 
served synchrotron radiation, and w, that of the 
electrons in the synchrotron. i indicates a phase 
difference of 7/2 between the two amplitudes, 
which means elliptical polarization with the limits 
of linear polarization at y=0 and circular polar- 
ization at very large angles where x = and K,, 
and K,, have the same asymptotic behavior. 

We have measured for both directions of polar- 
ization the variation of intensity with the angle 
with respect to the synchrotron plane, of the light 
emitted during the acceleration process at the 
Cornell synchrotron with peak energy 700 Mev. 
The apparatus used is shown in Fig. 1. The light 
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FIG. 1. Side view of light path through the apparatys, 


coming through the plane window in the vacuum 
chamber is focused by a lens with 7.8-inch focal 
length. Light emitted at the same angle y forms 
a horizontal line in the focal plane. There we 
have a horizontal slit, mounted on a vernier 
caliper, which enables fine vertical movement. 
Thus we select with the slit position certain 
angles of emission. To determine the direction 
of polarization we insert a polaroid filter as am- 
lyzer in front of the focusing lens. The light is 
detected with an RCA-6342 photomultiplier, the 
anode current being read directly on a micro- 
ammeter. At the end of the accelerating process 
the glass of the phototube fluoresces due to scat: 
tered electrons and their bremsstrahlung. The 
slit width and light filter are chosen to keep the 
signal large compared to this background. A 
yellow filter in front of the photomultiplier cuts 
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Photomultipiie, 








the spectrum off at about 4500 A; the upper limit 
is given by the response of the tube. The slit 
width was 4 mils or 0.5 milliradian. 

The result is shown in Fig. 2; the curve is the 
theoretical intensity with w =21c/(5000 A) and 
W, = 27 X10.8x10° sec™*. The difference between 
the theoretical curve and the experimental points 
is due to oscillation of the electrons, the band- 
width of the radiation, and the finite slit-width, 
all these effects being intrinsically more pro- 
nounced in the perpendicular component. The & 
perimental points are slightly asymmetrical, 
which is probably caused by a change of inclina- 
tion of the central orbit during the cycle. To in 
tegrate the total intensity in the two directions 0 
polarization, we measured the intensity with full 
width of the slit (1/8 inch). The result J,?/J,' 
= 3.45+ 0.08 is a little higher than the theoretical 
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FIG. 2. Experimental and theoretical intensities of 
the two polarization components of synchrotron radia- 
tion plotted as a function of the angle. 


value of 3. 

To check the phase relation between J, andJ,, 
aquarter-wavelength plate for the Hg 5461A line 
yas inserted between window and analyzer. We 
placed the slit at 11 milliradians and opened it 
48 milliradians. At this position the intensities 
vere J, = 19.5 units, J; =13 units. Inserting the 
quarter-wave plate and adjusting the analyzer at 
45° between the two directions, gave for the maxi- 
mum position 27 units, for the minimum 3. Thus 
the nearly circularly polarized light was con- 
verted into nearly linearly polarized light, which 


shows that the derived phase relation is correct. 

The measurement of the magnitude and the 
angular dependence of the polarization show good 
agreement with the theory in its general shape. 
Since our measurement integrates J,’ and J,” for 
one cycle and over all electrons, which might 
have an asymmetrical angular distribution 
changing with time, there are deviations in the 
details. The sensitivity to angular changes of the 
central orbit could enable one to make very fine 
beam studies. 

I am most grateful to Professor Robert R. Wil- 
son and the Staff of the Synchrotron for the hos- 
pitality extended to me. 
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PLASMA STABILITY AND BOUNDARY CONDITION 


F. C. Hoh 
Department of Electronics, The Royal Institute of Technology, Stockholm, Sweden 
(Received May 10, 1960) 


Recently, the validity of the classical theory 
of diffusion of charged particles across a mag- 
netic field has been extensively discussed.’ ~*° 
Inthe experiment made by Lehnert'»” agreement 
with the classical theory was found for magnetic 
fields weaker than a certain critical value Be 
beyond which this theory was observed not to be 
valid. The present note shows that this phenom - 
tion can be explained by an instability of the wall 
sheath. The experiments!» °-’,*° indicate that 
this instability causes the plasma to diffuse 
across the magnetic field much more rapidly 
than predicted by the classical theory. 

Bohm’s* criterion for a stable wall sheath 


. ( / ) , ( ) 


where v;’ is the ion velocity normal to the wall on 
entering the sheath, k is Boltzmann’s constant, 
T, is the electron temperature, and m; is the 
ion mass. In ordinary gas discharges’ °~’,™ it 
can be shown that this criterion remains approx- 
imately valid when the effect of a magnetic field 
(of the order of B,) is taken into account. 

Here it will be postulated that the plasma be- 
comes unstable when the boundary condition Eq. 
(1) is not satisfied because v;’ can be reduced 


559 






































VoLuME 4, NuMBER 11 





PHYSICAL REVIEW LETTERS 








by a magnetic field. 
For a cylindrical positive column in an axial 
magnetic field the radial ion velocity is *»° 


my k(T,+T,) (i (>) J, Qr/R) (2) 
.-@ B+B;) R J (ar/R)’ 





where 
B= (8, /P [1+ (BB ./b))), (3) 


and the corresponding equation when i is re- 
placed by e. Here \=2.4, R is the tube radius, 
ry is the distance from the axis, J, and J, are the 
Bessel functions of orders 1 and 0, e is the elec- 
tron charge, p, is the gas temperature at 0°C, 

T; is the ion temperature in the plasma, 8,; and 
8; are the electron and ion mobilities at unit 
pressure and at the corresponding temperatures, 
and B is the magnetic field strength. Combining 
Eqs. (1), (2), and (3) and putting r=r’, corre- 
sponding to the position of the sheath edge, we 
obtain the stability criterion: 


RS = -l/2 , e 2\V2 
B B. (84859) (g’p,/R by) , (4) 


where 


(5) 





“wT 47 r/ 
r-(7) ott; Bey 8jy %y0”'/®) 


T e BL4+8i, Jyiar’/RY 


Near the wall J,(Ar’/R) = 0.52. J,(ar’/R) =n,'/n, 
is the ratio of the electron density at r=r’ to 
that at the axis, and can be written as’ 


J 9 (r/R) = 





—_ , 2 , 
€_k TA, Ts, tee sad /R) a. 
B.+B, \ R n 


+s , 
ne 2 
0 


Here €, is the dielectric constant of vacuum, 
and n;’ is the ion density at r’. The position r’, 
where the plasma merges into the sheath, can 
be quite well defined by requiring” (n;’ -,’)/ng’ 
= 0.01-— 0.05, as B. is rather insensitive to this 
value [see Eqs. (4), (5), and (6)]. 

Earlier measurements’ *' of B, are plotted 
against (p,/R)” in Fig. 1. Under these discharge 
conditions we have on the average T, /V5= 1200 
°*K/volt (V; is the ionization potential), 7; ~400°K, 
and J,(Ar’/R) = 0.013 [Eq. (6)]. With the average 
mobility data, ** B,, is calculated from Eqs. (4) 
and (5) for helium (full curves in Fig. 1). Sim- 
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and Levine (reference 5) in the afterglow of a micro- 
wave discharge, compared with the present theory. 
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FIG. 1. Earlier measurements (references 7 and 1]) 
of B. as a function of (p)/R)”*. The curves are calcu- 
lated from Eqs. (4) and (5). 


ilar results are also obtained for discharges with 
krypton, ”»** argon, * >" and neon.® The present 
theory may also explain the critical magnetic 
field observed by Bostick and Levine’ in the 
afterglow of a microwave discharge (Fig. 2). 
Their experiment indicates that the onset of 
instability is due to wall conditions and is not 
associated with an axial electric field. 

The above discussions are valid for partially 
ionized gases. The corresponding situation for 
a fully ionized plasma remains to be investigated 
However, the initial stage of the discharge in 
the B-3 Stellarator*>® is not very far from that 
of the present experiments, ’»™ but with magnetic 
fields probably too strong to sustain a stable 
wall sheath. If the plasma is kept from the wall, 
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e.g., by means of apertures inserted in the dis- 
charge tube, the wall conditions are modified. 
This also seems to be supported by experiments.*® 
Further, the present theory seems to connect the 
classical diffusion theory’’* and the ion wave 
instability theory* in a natural way. 

An important consequence of the present dis- 
cussions is that a reversal of the radial electric 
field in a plasma by means of a magnetic field 
seems to be impossible if the boundary condition 
Eq. (1) has to be satisfied. Finally, for large 
values of Rp,, Eq. (4) shows that the sheath is 
unstable even in the absence of a magnetic field. 
This may be of some interest in the study of 
high-pressure glow discharges.** 

The author wishes to express his sincere thanks 
to Professor H. Alfvén, Dr. B. Lehnert, Dr. S. 
Lundquist, and Dr. E. Astrém for many valuable 
discussions. 
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ZEEMAN EFFECT IN THE RECOILLESS y-RAY RESONANCE OF Znf 


P. P. Craig, D. E. Nagle, and D. R. F. Cochran 
Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico 
(Received May 4, 1960) 


Recoil-free resonance emission and absorption 
(Mossbauer effect)’ of the 93-kev y ray in Zn” 
gives rise to the most precise energy definition 
thus far reported.? Despite the numerous diffi- 
culties which beset the experimenter searching 
for the resonance,°® a small but definite Moss- 
bauer effect has been found.? This Letter reports 
measurements on the influence of the nuclear 
Zeeman effect and other perturbing factors upon 
the Mossbauer effect in Zn™ embedded in an en- 
tiched ZnO absorber lattice. 

The relatively high energy of this y ray makes 
itnecessary to embed the source and the ab- 
sorber atoms in rigid crystalline lattices, and 
to perform the experiments at low temperatures. 
The first of these requirements was met by using 
ZnO for both the source and the absorber lattices. 
The second requirement was more than satisfied 
by using temperatures below the helium lambda 
transition (2.175°K). In addition to these basic 
requirements, several experimental difficulties 


are consequences of the extreme narrowness of 
the line (4.84x10-" ev). It is accordingly neces- 
sary to take into account the effect of various 
perturbing influences. Here we list the more 
important of these perturbations, and in the fol- 
lowing paragraphs indicate how they enter into 
the design of the experiment. 

One class of shifts arises from the change of 
nuclear.mass upon y-ray emission or absorption, 
with a resultant change in the phonon spectrum 
of the lattice. The change in energy of the emit- 
ted y ray is given by* 


AE = -(E/Mc*){T), (1) 


where E is the y-ray energy, M is the mass of 
the emitting nucleus, and 7’) is the expectation 
value for the kinetic energy per atom of the lat- 
tice. If any parameter x should differ between 
source and absorber lattice, the recoil-free 
peaks will occur at different energies in the 
emission spectrum and in the absorption spec- 
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trum. The shift is given by 
fe) 
5(AE) = = (AE) dx. (2) 


This expression predicts that a difference in 
Debye temperature between source and absorber 
of only about 1.3°K would cause a shift of one 
linewidth. Thus a change in the average iso- 
topic mass number of 2% would, through the 
mechanism of the Debye temperature (taking @ 


= 300°K), result in a shift of about four linewidths. 


Similarly, differences in chemical constitution® 
or lattice defects may be expected to produce 
significant shifts. We shall refer to a shift due 
to difference in isotopic mass between source 
and absorber as the isotopic mass effect. 

In order for the emitted gamma radiation to 
remain unshifted by the recoil of the emitting 
nuclei, it is essential that the recoil momentum 
be absorbed by entire crystallites. For all pre- 
viously observed resonances the minimum size 
of the crystallites was quite small. In the pre- 
sent situation this is no longer true. The recoil 
momentum must be taken up by at least 2 x 10° 
nuclei, so that the ZnO crystallites must be 
larger than 0.4 micron. Since commercial ZnO 
would normally possess average grain sizes 
smaller than one micron,® care must be exer- 
cised to ensure a grain size much larger than 
the above nominal requirement. A sintering pro- 
cess assured that this requirement was met for 
the source. A sample of the enriched absorber 
was studied under an oil immersion microscope. 
No grains smaller than 0.5 micron were ob- 
served, and the majority of the grains were in 
the range of 1 to 2 microns. 

Mechanical vibrations of only 10°° cm/sec 
would produce a Doppler broadening of about one 
linewidth. Since the vibration level in our build- 
ing was nearly 5x10~* cm/sec, the helium cryo- 
stat was shock- mounted and the pumping lines 
carefully decoupled. The source and the ab- 
sorber were clamped rigidly in a single package, 
which was suspended by threads in the helium 
bath. The helium bath was pumped below the 
lambda transition and the nitrogen radiation 
shield was frozen to prevent vibrations from 
boiling liquids and also to place source and ab- 
sorber in an isothermal bath. 

In place of the velocity drive usually used in 
this type of experiment, the resonance was 
shifted by means of the nuclear Zeeman effect. 
The magnetic moment of the Zn™ nucleus inter- 
acting with an applied magnetic field splits the 
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ground state into six components and the excited 
state into four components. Selection rules per- 
mit twelve component gamma transitions, as is 
shown in Fig. 1(a). Here /,=5/2, J, =3/2, and 
the ground-state gyromagnetic ratio g= 0.35.’ 
For the excited state no g value has been meas- 
ured; we adopt arbitrarily the value -1. (The 
negative sign is predicted by the shell model.) 
The source is shielded from the magnetic field 
and hence the emission spectrum is not split. 
For simplicity, the effect of quadrupole interac- 
tions has been omitted. Assuming a shift 5E be- 
tween source and absorber, resonance can occur 
for six values of magnetic field. Each component 
has a width which is compounded of the natural 
linewidth and the widths due to residual mechan- 
ical vibrations, quadrupole broadening, etc. 
Even if this width is the same for all components, 
the present method of studying the resonances as 
a function of an external magnetic field will 
cause the apparent or magnetic width of each re- 
sonance peak to be proportional to the field re- 
quired to establish the resonance. Thus reso- 
nances occurring at low values of the applied 
field will appear narrower than those established 
at high fields, and the observed spectrum will 
appear distorted. 

The magnetic field was produced by a small 
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FIG. 1. (a) Zeeman splittings of the gamma ray vs 
an applied magnetic field. The gyromagnetic ratio of 
the ground state is 0.35, and a value of -1.0 has been 
arbitrarily chosen for the excited state. (b) Resonance 
pattern expected from the splittings in (a) in the pre- 
sence of an energy shift 6E between source and ab- 
sorber. In actuality, line broadening would be ex- 
pected to merge closely spaced lines into single peaks. 
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solenoid surrounding the absorber. Fields of up 
to 700 gauss could be applied to the absorber, 
while the source was entirely shielded (at oper- 
ating temperatures) by superconducting lead 

foil. Stray fields at the absorber were measured 
to be less than 2 gauss. Fields:from the sole- 
noid were prevented from reaching the photo- 
multiplier by a superconducting ring placed in 
the cryostat bottom through which the gamma 
rays emerged. Changes in counting rate due to 
changing fields at the photomultiplier were found, 
even without the superconducting shield, to be 
less than + 0.02 %. 

The absorber used in all experiments was 
1.231 g of enriched ZnO (92.4% Zn®™ enrichment) 
pressed with 0.036 g of polyethylene glycol bind- 
er into a button 1.11 cm in diameter. Sources 
were prepared in the Los Alamos cyclotron by 
the reaction Zn®(d,n)Ga® on normal sintered 
ZnO. Other reaction products were accounted 
for by a background correction. After bombard- 
ment, the sources were annealed for about one 
hour in air at about 1000°C in order to anneal 
radiation damage and to assure that the Ga®”’ 
atoms were correctly placed in the ZnO lattice. 
Such a procedure is essential, and lower tem- 
peratures or shorter times yielded erratic re- 
sults, or no resonance whatever. 

Measurements were made of the transmission 
of the absorber vs the applied magnetic field. 
Because of the extremely small change in count- 
ing rate and the large number of counts required 
(typically 10° per point), automation of the count- 
ing system and extreme system stability were 
imperative. The 93-kev gamma rays passed 
through thin windows in the cryostat bottom, and 
were detected by a Nal(T1) scintillation crystal. 
The pulses were amplified and analyzed by a 
single-channel analyzer followed by a scaler 
equipped with a digital recorder. In order to 
eliminate the effect of source decay, all meas- 
urements were taken relative to the counting 
rate when one applied a magnetic field sufficient 
to destroy the resonance entirely. The magnetic 
field sequencing and the scaler were controlled 
by an automatic programmer. Timing was con- 
trolled by a thermostatted quartz crystal fre- 
quency standard with a stated accuracy of 3 parts 
in 10’ per week. The programming was such as 
to make the measurements insensitive to linear 
drifts. 

Instrument checks were performed by cycling 
the magnetic field (a) at room temperature, (b) 
at low temperature using a nonresonant gamma 


ray, and (c) at low temperature using as an ab- 
sorber a zinc-containing material (gahnite) in 
which no resonance exists.? In each case the 
change in counting rate upon application of the 
magnetic field was found to be zero to within 

+ 0.02%. 

Figure 2 shows the resonance curve obtained. 
The points include corrections for background 
and unresolved nonresonant gamma rays of 
typically 30%. The central features of the data 
are the remarkably large total area under the 
resonance, and the structure and total breadth 
present. The total area is sufficiently large that, 
were there no line broadening, one would esti- 
mate the resonance absorption to be several per- 
cent. Such a value implies a Debye temperature 
for ZnO of about 300°K, which is consistent with 
the results obtained from specific heats. The 
major structural feature of the curve lies in the 
fact that the maximum resonance does not occur 
at zero field, but is shifted to about 10 gauss. 
This result indicates the presence of shifts such 
as the isotopic mass effect mentioned above. 
However, the displacement of the maximum (to 
10 gauss) is smaller than one would predict from 
this effect alone unless an unreasonably large 
value is assumed for the excited state g value. A 
plausible explanation is that chemical shifts® are 
superimposed upon the isotopic effect. The pre- 
sence of such shifts is implied by the small re- 
sonance found at zero field using an unenriched 
absorber.” Since in that case no isotope effects 
were present, only chemical shifts or quadrupole 
broadening can explain the small resonance ob- 
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FIG. 2. Experimental resonance absorption pattern 
for the 93-kev line of Zn*’ vs the magnetic field applied 
to an enriched ZnO absorber. The source is in zero 
magnetic field. I'/u, is the natural level width (in 
ergs) divided by a nuclear magneton (in ergs per 
gauss). 
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served. In the present experiments two peaks 
are found separated by only about 40 gauss so 
that it seems unlikely that the broadening is ex- 
treme. The tail of the resonance extends to 
about 400 gauss with no additional structure. 
This may be related to the characteristic of the 
magnetic method that even unbroadened lines 
possess apparent magnetic widths which are pro- 
portional to the applied magnetic field. 

Although the interpretation is admittedly in- 
complete, the extreme sharpness of the reso- 
nance is apparent. In further study, involving 
the development of a Doppler shift drive, we 
hope to measure a number of the energy shifts 
and level splittings mentioned in previous para- 
graphs. 

We wish to thank S. D. Stoddard and R. E. 
Cowan for preparation of the ZnO source buttons 
and for compacting the enriched ZnO absorber. 
The generous cooperation of the cyclotron group 
is gratefully acknowledged. W. E. Keller and 


J. G. Dash each contributed a number of ideas to 
the experiment. 
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OPTICAL AND MICROWAVE-OPTICAL EXPERIMENTS IN RUBY 


T. H. Maiman 
Hughes Research Laboratories, Malibu, California 
(Received April 22, 1960) 


Several recent papers’ * have reported optical 
and microwave-optical measurements in ruby 
(Cr*** in Al,O,). We wish to report here some 
new experiments concerning the fluorescent re- 


laxation processes in this crystal. Reported here 


also are the first observations of ground-state 
population changes in ruby due to optical excita- 
tion and the detection of optical absorption be- 
tween two excited states in this crystal. 

The predominant processes which ensue in a 
fluorescent material when it is irradiated at an 
appropriate wavelength are shown in Fig. 1. W,, 
is the induced transition probability per unit time 
due to an exciting radiation and the S,,, are de- 
cay rates which include both radiative and non- 
radiative processes. In this crystal S,, is easily 
obtained from the decay rate of the fluorescent 
level (EZ) after an exciting source is turned off. 
The lifetime for this process is about 5 msec. 
Varsanyi, Wood, and Schawlow’ have further de- 
monstrated that this lifetime is almost entirely 


due to spontaneous emission, i.e., S,, is approxi- 


mately the Einstein A coefficient A,,. 
An approximate value for the rate S,, was ob- 
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tained in the following way. A crystal of ruby 
was irradiated with 5600A radiation causing ab- 
sorption into the lower band (*A,—‘*F,). The sam- 
ple used was a one-centimeter cube cut from a 
boule of standard pink ruby supplied by the Linde 
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FIG. 1. Pertinent features of the energy level dia- 
gram in ruby. 
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Company, with a concentration of approximately 
0.05 weight percent of Cr,O, to Al,O,. Two com- 
ponents of radiation re-emitted from the crystal 
were observed in a direction perpendicular to 
the exciting beam: that due to re-emission of 
the incoming radiation (spontaneous decay from 
‘F,) and fluorescence (spontaneous decay from 
FE), The intensity of the first component is pro- 
portional to hv,,N,A,,, where A,, is the A coeffi- 
cient for *F,~*A, and is calculated from meas- 
urements of absorption coefficient and linewidth 
for this transition (A,,~3x10°/sec). Similarly 
the fluorescent intensity is proportional to 
hv,,N,A,,- By a measurement of the ratio of 
these two components and the use of an auxiliary 
condition applicable to steady-state conditions 
N,S,,=N,S5, and also the use of the approxima- 
tion S,,=A,,, we find S,,~2x10"/sec. 

A measurement of fluorescent quantum effi- 
ciency, i.e., the number of fluorescent quanta 
emitted compared to the number absorbed by the 
crystal from the exciting beam, yielded a value 
near unity. This result reconfirms the evidence 
that the life of level 2 is near radiative and also 
implies that S,,>>S,,. The experiment was not 
accurate enough to yield a precise value but does 
indicate that the nonradiative process (S,, - A,,) 
<4x10°/sec. 

Calculations utilizing the previous results in- 
dicated that population changes in the ground 
state of ruby due to optical excitation would be 
easily observed. This conclusion was verified in 
the following experiments. A ruby crystal was 
mounted between parallel silvered plates to form 
a microwave cavity resonant at the ground-state 
zero-field splitting (11.3 kMc/sec). About half 
the cavity losses were due to magnetic absorption 
as evidenced by an increase in cavity Q when a 
small magnet was brought near the ruby. The 
reflection coefficient of the cavity was monitored 
on an oscilloscope while a short pulse (200 psec) 
of light from a flash tube irradiated the crystal. 
The magnitude of the microwave magnetic absorp- 
tion was observed to decrease abruptly and then 
return to equilibrium with a time constant of 
about 5 msec (see Fig. 2). We attribute this ef- 
fect to temporary depletion of this ground state 
population with subsequent decay back from the 
fluorescent level. The experiment was performed 
at room temperature where the thermal relaxa- 
tion times in the ground state of ruby are the 
order of a microsecond; in the time scale of the 
experiment, therefore, Boltzmann equilibrium 
in these levels is maintained. 








FIG. 2. The upper trace was obtained from a photo- 
multiplier monitoring the R-line fluorescence (7E—‘A,) 
in ruby. The lower trace is a recording of the reflec- 
tion coefficient of a microwave ruby cavity during the 
same time interval. The time scale is 2 milliseconds 
per division. Both recordings followed an intense light 
pulse which excited the crystal via ‘A,-‘F,. 


A repetition of the above experiment at liquid 
helium temperatures is being planned. At this 
temperature we would expect to be able to ob- 
serve directly any preferential depopulation of 
the ground sublevels due to polarized light and 
also any preferential repopulation of these levels 
since the thermal relaxations times would then 
be 30-100 msec. 

To verify further the depletion of ground-state 
population observed in the previous experiment 
an independent measurement was made. A beam 
of monochromatic light of wavelength 4100 A was 
transmitted through a ruby crystal and partially 
absorbed due to the transition “A,—*F,. When the 
intense pulse of radiation at 5600 A was turned 
on, the 4100A radiation passing through the crys- 
tal abruptly increased and subsequently decayed 
in about 5 msec just as the microwave signal in 
the previous experiment. This result was ex- 
pected since the temporary reduction in ground- 
state population caused the crystal to become 
more transparent to the 4100A radiation until the 
fluorescent level decayed to normal. In both ex- 
periments a population change of about 3% was 
estimated. 

An unexpected result was observed when the 
probe wavelength was changed from 4100 A to 
3600 A. In this case a decrease in light intensity 
emerging from the crystal was observed. This 
implies that the crystal became more absorbing 
even though the ground-state population was de- 
creased. We can explain this last effect when it 
is realized that 3600A radiation can cause transi- 
tions from the fluorescent level (?E) to a high- 
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lying charge transfer band (not shown in the fig- 
ure). Consequently, we conclude that we were 
observing transitions between two excited optical 
states. The fact that the abruptly increased 
3600A absorption also decayed with a 5-msec 
time constant is consistent with and strengthens 
the above conclusion. 

The author gratefully acknowledges helpful dis- 
cussions with G. Birnbaum, R. W. Hellwarth, 
L. C. Levitt, and R. A. Satten and is particularly 









indebted to I. J. D’Haenens for technical assist- 
ance in performing the experiments. 
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Rosenberg, Spruch, and O’ Malley’ (RSO) have 
recently calculated rigorous upper bounds, 
a, <6.23 and a;<1.92, for the singlet and triplet 
scattering lengths in the scattering of electrons 
from atomic hydrogen,” under the assumption 
that the H ion has only one bound state, with 
singlet spin. The fact that these upper bounds 
are 15% below the results of most of the elabo- 
rate calculations thus far performed’ indicates 
that we are far from an adequate quantitative 
theory of this three-body problem. 

Those calculations which yield scattering 
lengths close to these upper bounds make use 
of the adiabatic hypothesis, according to which 
the distortion of the atom follows the instantan- 
eous position of the scattered electrons.*»> Ob- 
jections to the validity of the adiabatic hypothesis 
have recently been raised.*° We shall outline in 
this note a rigorous (and hence nonadiabatic) 
series for the phase shift. The proposed series 
serves two purposes: First, because it is rapid- 
ly convergent it provides a method for calcula- 
ting phase shifts reliably. This is so because 
methods can be given for accurately approxima- 
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ting the higher order terms. Second, the suc- 
cessive terms of the projected series can be 
put in a one-to-one correspondence with an adi- 
abatic series which has been derived for dealing 
with this problem.**” It can be demonstrated 
that the main contribution to the nonadiabatic 
terms comes from the functions which enter the 
adiabatic development. This provides a justifi- 
cation for the qualitative validity of the adiabatic 
theory and an explanation of the fact that the 
adiabatic calculations give scattering lengths in 
better accord with the RSO upper bounds. 

The method is an extension of a method used 
originally by Breit and more recently by Luke, 
Meyerott, and Clendenin® for calculating the 
energy of some excited states of 2-electron atoms 
and ions. One makes a decomposition of the zero 
angular momentum wave function ¥(7,7,4,,), 


¥(7) 79949) 








_—— > (21+1)76 (7,75)P (cos®,,), (1) 
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and reduces the Schrédinger equation to the in- 
finite set of coupled equations (for 7, > 7.) 
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The asymptotic form of © for large 7, in the 
scattering problem is such that all the $7 for 
1>0 vanish. For #, we have 


lim $5 (747%») =sin(kr, +O)R, (r 3 (3) 


§ is the required phase shift, and R,,(7,) is 7, 
time the radial ground-state wave function of the 
hydrogen atom. 

Now consider the equation obtained by neg- 
lecting the right-hand side of the /=0 equation: 


e” + ¥ onagel }, 7,72) =0, 7,272. (4) 
or,” er,” % 0 "as ’ 7 


We seek a solution of this equation in which 
6,” has the same asymptotic form as (3), but 
with a phase shift 5, in genera: different from 6. 
The main formula of this approach, which may 
readily be derived from (2) and (4), relates 6 
and 6,: 








ksin(6 - 5,) 
[2] ” eo v; Y. l 
= -F) 2(21 +1) ay | ar, dr,6,°2_# . (5) 
4 piel 
0 0 


The demonstration that the successive terms 
on the right-hand side of (5) become rapidly 
small depends on showing that the major contri- 
bution of each term comes from the region 
Y,“yr,. That this is so follows from the asymp- 
totic form of #) in that region, 


-2sin(kr, + 5) = F rn) 
1 > 


im ¢,° 1/2 141° mn 
Y 


n= ? (9744) 





I>0. (6) 


This satisfies (2) when one neglects all terms 
proportional to r,~’, where n >1+2. It also 
indicates that @7, for small 7,, vanishes as an 
inverse power of 7,, as opposed to its exponen- 
tially decreasing behavior for rv, and 7, both 

large. [The right-hand side of (6) is essentially 
the /th term of an adiabatic series which was pro- 
posed in references 3 and 4 for dealing with this 
problem.] In the region where r, and 7, are both 
small, } j is itself small by virtue of the centri- 
fugal potential term in (2). The diminishing of 
the terms in (5) as a function of J can be explicit- 
ly shown by using (3) and (6) for 4, and #). 

The physical meaning of (4) corresponds to 
electron 1 coming in from infinity and seeing no 
force at all until it is actually inside electron 2, 
at which point it sees the complete nuclear 















Table I. Approximate 6, (in radians) based on two 
terms in the expansion of ,™. 





A 0.05 0.1 0.3 0.5 0.75 





5p (singlet) 2.76 2.41 1.48 0.92 0.41 
59 (triplet) 3.03 2.93 2.49 2.09 1.71 





charge. Or to put it another way, when electron 
1 gets inside electron 2, the latter becomes the 
scattered particle. Equation (4) is only valid 

for r,>7,, and the symmetry of 4, is taken 
into account by demanding that 4, or its nor- 
mal derivative (corresponding to triplet or sing- 
let scattering) be zero along the line 7,=7,. The 
totality of separable solutions of (4) can readily 
be written down. In Table I we have given 6, 
obtained from including two such terms so that 
their sum “best” satisfies the boundary condi- 
tion along 7, =7,.° 

The scattering lengths obtained from the k=0.05 
entries in the table are a, =7.6 and a;=2.2. These 
are above the Spruch-Rosenberg limits, and they 
are similar to the results of calculations using 
correctly symmetrized wave functions that do 
not naturally describe a long-range polarization 
of the atom, the exchange approximation for 
example.'® Our main point is that the inclusion 
of the first term of (5) will bring the scattering 
lengths near or within the Spruch-Rosenberg 
limits. This has already been indirectly tested 
by the method of polarized orbitals* in which the 
exchange wave function is modified by a polari- 
zation term, which is essentially the asymp- 
totic form (6) of #,, and which reduces the scat- 
tering lengths below the RSO upper bounds. 

One is, of course, interested in the phase 
shifts as well as the scattering lengths (and not 
only in upper bounds for these quantities). The 
relevance of the higher terms in (5) is that they 
allow one to calculate the phase shifts in a con- 
vergent manner. To all these calculations, the 
solution @, of (4) is basic. For if that is known 
accurately, the contribution of the other terms 
can be reasonably approximated from the asymp- 
totic form (6) of p modified for smaller values 
of y, and ry, with parameters which can be de- 
termined by additional sum rules relating the 
$, to ,. 









“Most of this work was performed while the author 
was at the National Bureau of Standards. 

'L. Rosenberg, L. Spruch, and T. O’Malley, Phys. 
Rev. (to be published). See also L. Spruch and L. 
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METHOD FOR DETERMINING THE ORBITAL ANGULAR MOMENTUM IN K -d CAPTURE* 


Donald H. Miller 
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(Received May 9, 1960) 


In a recent communication, Day, Snow, and 
Sucher predicted that K” mesons when stopped in 
liquid hydrogen or deuterium would be captured 
from high-lying S atomic orbitals of the K -p or 
K’-d system.' They found that the Stark- effect 
mixing of angular-momentum states induced by 
the strong electric fields experienced by the neu- 
tral K-p or K -d system as it passed within the 
Bohr orbitals of other atoms of the liquid allowed 
rapid capture through the S-wave K-nucleon in- 
teraction. Their prediction has been used to draw 
important conclusions about the properties of 
hyperons’ and the nature of the low-energy, S- 
wave K-nucleon interaction.” In addition, esti- 
mates indicate that Stark-effect mixing deter- 
mines the course of the antiproton annihilation 
reaction at rest.* Because of the wide applica- 
bility of the prediction, an experimental test is 
highly desirable. The obvious experiment of look- 
ing for K-shell x rays when K” mesons are 
stopped in liquid H, or D, is technically difficult 
and probably will not be done in the near future. 

It may be noted that several of the K -d cap- 
ture reactions are sensitive to the orbital angular 
momentum of the K -d system. This is implied 
for example in the calculations of Pais and Trei- 
man* and Day and Snow’ on the = m-hyperfrag- 
ment production rate. Their work shows that if 
the bound =n system existed and its character- 
istics were known, the production rate for 
stopped K” mesons would provide a clear deter- 
mination of the orbital from which capture oc- 
curs. It is the purpose of this note to point out 
that the K -d orbital angular momentum deter- 
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mines the rates for another class of reactions 
that is more accessible experimentally, the non- 
mesonic capture reactions in deuterium, 


K° +d-= "+p, (1a) 
K +d~°+n, (1b) 
K~ +d~A+n, (1c) 


where the rate for reaction (la) is twice that for 
(1b) by charge independence. We shall confine 
our attention to reaction (1a) since it is easily 
recognized in a deuterium-filled bubble chamber 
and is known to occur at a rate of ~0.7% when 
stopped K mesons are captured.® For K™ cap- 
ture at rest, the center-of-mass momentum of 
the 5"p system is 511 Mev/c, and production is 
inhibited unless the K’ is absorbed while the two 
nucleons are within a distance ~0.4x107** cm. 
The small size of this effective interaction vol- 
ume implies that the rate for (1a) will be propor- 
tional to | d, g(r =0)!? or to |Vo,, p(r = 0) |? depend: 
ing upon the atomic orbital from which capture 
occurs. On the other hand, mesonic absorption 
can occur over the entire volume of the deuteron 
and will proceed through the S-wave K -nucleon 
interaction for capture from either S or P atomic 
orbitals.° Therefore, if capture occurs predom- 
inantly from S orbitals as predicted by Day, 
Snow, and Sucher, the fraction of nonmesonic ab- 
sorptions observed at rest and in flight must be a 
continuous and slowly varying function of K' mo- 
mentum. 

Making the explicit assumption that the beha- 
vior of the nonmesonic transition amplitude is 
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determined by the centrifugal barrier in the ini- 
tial state, and that only S- and P-wave absorp- 
tion is important at low momentum, we have’ 


ae (q) =(A + Bq’)pen. (2) 


Here g is the center-of-mass momentum of the 
initial K -d system in units of up 4c, where ung 
is the K -d reduced mass and c the velocity of 
light; p is the final center-of-mass momentum 
measured in units of its threshold value, p,=511 
Mev/c; and N is the atomic density of liquid deu- 
terium. The quantities A and B are to be deter- 
mined from measurement of the in-flight produc- 
tion cross sections. The capture rate from 7S 
ornP orbitals may then be estimated by using 
nS 


= - 2 
t. =cAlo, o(r=0) | ' (3a) 


P 
he =cB(h/uy go) V9, p(r=0) |?. (3b) 


Since the absolute transition rate cannot be 
measured, Eqs. (3a) and (3b) must be compared 
with the mesonic absorption rate. To obtain the 
latter, we assume that the absorption rate in 
deuterium is essentially the sum of the absorp- 
tion rates for the free neutron and proton.* The 
mesonic absorption rate from an ml orbital is 
then given by 


h b.+3b 
nl “aa 9 
c. =41c (—)h 2 ] ( | 47) Ps (4) 





where xp is the reduced K -p mass, and the 
average over the nucleon density distribution is 
evaluated by using the Hulthén form for the deu- 
teron wave function. The quantities 5, and 6, are 
the imaginary parts of the T=0 and T=1 zero- 
energy scattering lengths evaluated by Dalitz and 
Tuan.* Dividing Eqs. (3a) and (3b) by Eq. (4), we 
dbtain for the fraction of nonmesonic absorptions 
from S or P atomic orbitals: 


R>(0) = 0.123 x107* A, (5a) 


RF (0) = 0.0241 x10" B. (5b) 


At present, only the fraction of nonmesonic ab- 
sorptions occurring at rest is known. Since 

tither Eq. (5a) or (5b) should essentially account 
for the observed value, R(0) = 0.007, A or B may 
te determined and corresponding lower limits 

‘or Tm (q) estimated from Eq. (2). To obtain 
wer limits for the fraction of nonmesonic ab- 








Table I. Lower limit for fraction of nonmesonic ab- 
sorptions expected in flight when K capture at rest 
occurs from either S or P atomic orbitals. 








K™~ laboratory S-orbital P-orbital 
momentum (Mev/C) capture capture 
50 0.007 0 
200 0.019 0.016 
300 0.017 0.032 





sorptions expected in flight, I,,,,(q) may be 
divided by the mesonic absorption rate as given 
by 


r (a) = a (82) N, (6) 


where 0, and o, are the absorption cross sections 
for the T=0 and T=1 states of the K -p system.° 
Typical results are listed in Table I for the two 
assumptions of S- or P-orbital capture at rest. 

The ease with which the experiment can yield a 
definitive result depends upon the extent to which 
A is not equal to B in Eq. (2). In the analogous 
reaction 1*+d— pb+p, Bis much greater than A 
because of the anomalously small S-wave, 17- 
nucleon interaction. In the present case, the S- 
wave, K -nucleon interaction is very large, and 
if A is much greater than B, a measurement of 
low statistical accuracy would indicate that or- 
dinary collisional and radiative de-excitation 
mechanisms were inadequate to account for the 
observed nonmesonic rate at rest. 





"This work was done under the auspices of the U. S. 
Atomic Energy Commission. 

'T, B. Day, G. A. Snow, and J. Sucher, Phys. Rev. 
Letters 3, 61 (1959). 

*R. H. Dalitz and S. F. Tuan, Ann. Phys. (to be 
published). 

3Bipin Desai, Phys. Rev. (to be published). 

‘A. Pais and S. B. Treiman, Phys. Rev. 107, 1396 
(1957). 

5T. B. Day and G. A. Snow, Phys. Rev. Letters 2, 
59 (1959). 

®Data of N. Horwitz, D. Miller, and J. Murray, re- 
ported by L. Alvarez at the Ninth Annual International 
Conference on High-Energy Physics, Kiev, Russia, 
1959 (unpublished). 

"The general properties for such a transition ampli- 
tude have been discussed by K. A. Brueckner and 
K. M. Watson [Phys. Rev. 86, 923 (1952)] and K. M. 


Watson [Phys. Rev. 88, 1163 (1952)] in connection with 
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the reaction p+p—n’ + d. 

8It may be expected that final-state interactions will 
redistribute the reaction products among the available 
final states but will not seriously affect the absolute 
transition rate. T. B. Day, G. A. Snow, and J. Su- 
cher, in University of Maryland Technical Report 167, 
March, 1960 (unpublished), have estimated the effects 
of multiple scattering on the absorption rate at 200 
Mev/c K~ laboratory momentum and found them to be 


small. 

%At zero and 200 Mev/c, the necessary parameters 
for the T=0 and T=1 absorption channels were taken 
from reference 2 for the experimentally preferred 
(a+) and (5+) solutions. At 300 Mev/c, we used the 
dy and o, absorption cross sections of P. Norden, R. 
Tripp, A. Rosenfeld, and F. Solmitz reported by L. 
Alvarez at the Ninth International Conference on High- 
Energy Physics, Kiev, Russia, 1959 (unpublished). 





FORMATION OF ».-MESONIC MOLECULES IN H-D MIXTURES* 


J. G. Fetkovich, T. H. Fields, G. B. Yodh, and M. Derrick 
Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
(Received May 9, 1960) 


Since the discovery by Alvarez et al.’ of the 
catalysis of nuclear reactions by negative muons, 
first discussed by Frank,’ there have been nu- 
merous estimates*° of the rates of the many pro- 
cesses involved. The results of these calcula- 
tions are disparate, even in some rather careful 
recent work.®,” Apart from the intrinsic interest 
in mesonic atomic and molecular systems, an 
additional stimulus to their study is the necessity 
of understanding these phenomena for the inter- 
pretation of future measurements of the rate of 
absorption of negative muons in hydrogen. 

Herein we report new experimental data, which 
allow a determination of some of the reaction 
rates involved in the catalysis processes. 

In the experimental work done previously, 
mesons were stopped in hydrogen containing 
small quantities of deuterium. The chain of 
events believed to lead to nuclear catalysis in 
that case is depicted in Fig. 1. Within this 
scheme, the number of nuclear reactions 


p+d—He*®+Q (1) 


(where Q signifies a y ray or a “regenerated” p~) 
per muon is given by 


Nia7'p r pd! “nHt D*pa 


x [rv Ap’ Onn*> 


+ Ao) 
tra @) 


where Spd is the fraction of (pu ~d) molecules 
which undergo nuclear reactions, and Cy, is the 
deuterium concentration (the hydrogen concentra- 
tion is 1-C))=1 for these cases). In deriving 

Eq. (2) from Fig. 1 the branch labeled CpApp is 
ignored, as is that labeled CpA,,--" In addition, 
recycling of the muon is ignored since it is ex- 
pected* to occur only for regenerations. This 
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makes an error of ~3% in N. 

The data so far available consist of the (bubble 
chamber) measurement by the Berkeley group’® 
of aNd [a is the fraction of all p-d nuclear re- 
actions leading to a regenerated muon: a =k/(k+1), 
where & is the internal conversion coefficient], 
and the experiment of Ashmore et al.*® measuring 
Npd at Cp=0.018. The bubble chamber data show 
that aN» is saturated at C,)=0.01 where aNpgq 
= 0.024+0.002. Furthermore, if one inserts the 
Berkeley data into Eq. (2), one obtains the rela- 
tion*» 





d,)- (3) 


= (2.0+ 0.7) x100..7+ 0 


r pd 
Be 
FAST 


(2p) —A0_. wecay) 


Copa 
(ud) —A0_» 


foo\on, 


(dpa) I pod 


CHAHH ‘nee 


(Pep) toure 


Xo Xo 





A. 
REACTION (!) 


NUCLEAR 
REACTION (1) 


FIG. 1. Chain of events believed to lead to catalysis 
at low deuterium concentration. Note the scheme for 
labeling and ordering the reaction rate subscripts. 
The reaction rates are defined as those for pure hydro- 
gen and deuterium, at bubble chamber conditions: 
pyH= 0.58 g/cm’, and pp=0.133 g/cm’, The double 
lines show those processes which are important at 
saturation. We neglect recycling, as mentioned in the 
text. 
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The experiment of Ashmore et al. indicates" that 
at saturation, 


7.4" 0.344 0.06, 


0.55 x10 sec™*<A,,< 2 x10" sec™’, 


0.30<f,5< 0.66. (4) 
In the present experiment, we stopped negative 
muons in a bubble chamber containing 94.9 atomic 
% deuterium and 5.1 atomic % hydrogen. The ex- 

pected chain of events in this case is shown in 
Fig. 2. Again the fraction /, d of (pu~d) molecules 
undergo reaction (1), while the fraction fgg of 
(du-d) molecules undergo one of the reactions 


d+d~He*+n, (5) 
d+d-~—H*+p. (6) 


Assuming complete recycling of the muon in re- 
actions (5) and (6),* the number, Nig, of nuclear 
reactions of type (5) or (6) per muon is 


aa ~faa© ppp! py" Saga pp* wat *o! (7) 





p< 


FAST 


id) ———— >» 


he \~ 


at 9,- d) oo — 2» 


v 
NUCLEAR 
REACTION (5),(6) 


NUCLEAR 
REACTION (1) 











FIG. 2. Chain of events believed to lead to catalysis 
at low hydrogen concentration. Recycling from re- 


action (1) is neglected, while recycling in (5) and (6), 
expected to be 92% (see reference 4), is assumed to 
be complete. 






while the number of regenerations per muon is 


ON a= Spa ww / [op Sa@pp* Cw Dat *o! 
(8) 


In this experiment we found that the number of 
regenerations per muon is nearly equal to the 
number found by Alvarez et al. at saturation, 
indicating that App is quite large. The events 
found were (including all necessary corrections) 
96+ 16 cases of reaction (6) [we cannot detect re- 
action (5)] and 33.84 5.6 regenerations in pictures 
which contained a total of 2201+ 50 y -e decays. 
Reaction (6) and the regenerations were identified 
by the characteristic ranges of the proton and the 
muon, respectively. In addition to these events, 
62+ 8 mesons were absorbed in the chamber. We 
concluded from a counter investigation of the 
beam as well as curvature measurements in the 
chamber that these were not pions, and therefore 
must be attributed to muon capture on impurities 
in the chamber. 

From the ratio of the number of cases of reac- 
tion (6) to the number of stopping muons we see 
that {7 7 0.1 (it is expected* that fgg=1), where 
we have assumed, as we will throughout, a 
branching ratio of 1:1 between reactions (5) and 
(6)..* We could in principle, by inserting our 
data into Eq. (8) and comparing this with Eq. (2), 
using the Berkeley saturation data, determine 
ApH: However, due to uncertainty in the Z of the 
impurity, we cannot definitely exclude ApH = * 
Nevertheless we can obtain a lower limit to ADH 
by assuming that all impurities have high Z, and 
by assuming that f7,=1. We thus obtain 


>1x10" sec™* (90% statistical probability), 
(9) 
OF "0-084 s 0.008. (10) 


‘DH 


Further, by dividing Eq. (7) by Eq. (8) we obtain 


(a. /A,.,) = (8.54 3) x 1073, (11) 


Sad DD DH 


which is independent of the Z of the chamber con- 
taminant. Equations (11) and (9) yield 


Sad pp? 8:5 *1% sec™*. (12) 


Inclusion of Ashmore’s results, Eq. (4), indi- 
cates that A,,,, =2 10" sec"! and fq = 0.34, and 
thus that fggApp =1.710° sec”. For comparison, 
a recent calculation® yields the values ADH =2.5 
x10® sec” and App* 5.9x10* sec™}. 
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The large discrepancy between the calculated 
and experimental values is surprising. If the ex- 
periments are correct, something is wrong with 
the calculations of meson- molecular formation 
rates, or the mechanisms (Figs. 1 and 2) sup- 
posed to lead to nuclear catalysis are incorrect. 

The analysis given here indicates that in liquid 
hydrogen containing >1% deuterium, all absorp- 
tions of »” by protons will be from a (pu d) 
molecule. For liquid hydrogen containing <1% 
deuterium, it is not at present possible to know 
the molecular or atomic states from which muons 
are captured. However, even for absorptions 
from (pu-d) molecules, capture and absorption 
of the muon by the He® reaction product will seri- 
ously complicate the interpretation of experi- 
mental data. 

Work is presently in progress on an extension 
of this experiment, with a smaller concentration 
of hydrogen, which should allow a more precise 
determination of » and X D’ 

We would like to tfank Gate Pewitt, John Deahl, 
and Leo Fatur for assistance in obtaining the 
data, and Anna Marie Bankowski and Donna Sta- 
sak for scanning the film. 
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HIGH-ENERGY BREMSSTRAHLUNG FROM A SILICON SINGLE CRYSTAL 


G. Bologna, G. Diambrini, and G. P. Murtas 
Laboratori Nazionali di Frascati del Comitato Nazionale per le Ricerche Nucleari, Frascati, Italy 
(Received May 3, 1960) 


In a previous Letter’ we reported the results 
relative to electron pair production from a sili- 
con single crystal. In this Letter we give the re- 
sults relative to several measurements on brems- 
strahlung from a similar target. We also com- 
pare these results with the theoretical prediction 
by Uberall? and Schiff.° 

We used about the same experimental arrange- 
ment described in reference 1. The only differ- 
ences are that the single crystal is now mounted 
within the synchrotron chamber and that the 
spectrometer converter is an aluminum one; 
moreover we added another counting channel, for 
measuring simultaneously at two different photon 
energies. 

The silicon single crystal is in the form of a 
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half-circular plate 15 mm in diameter and 2.7 
x10°* radiation length in thickness; it is cut per- 
pendicular to the axis [111] within + 4 mrad, as 
determined by a Laue x-ray back-reflection meth- 
od.* A goniometric device allows the single crys- 
tal to be rotated both about a horizontal and a 
vertical axis; the precision in the measurement 
of the angles is + 0.5 mrad. 

First we successively measure the numbers 
N(@,k), N(@,k,) of symmetrical pairs per fixed 
number of monitor units (corresponding to 10*° 
equivalent quanta), as a function of the angle @ 
between the incident 1-Gev electron beam and 
the crystal axis, and for the central values ,?, 
of the photon energies. We subtract for delayed 
coincidences and for background as in the pre- 
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vious work.’ 

We choose k, = 900 Mev because, as follows 
from Uberall’s theory,” the bremsstrahlung inten- 
sity at such an energy has a negligible depend- 
ence on 6, so that N(6,k,) may be used as a nor- 
malization factor. The y-ray beam monitor units 
cannot be used, because the shape of the spec- 
trum depends on 6. N(6,k,) ranges from ~ 8000 
to ~ 12000 counts/monitor unit. 

Due to the multiple traversals of the electrons 
through the crystalline target, the y-ray beam 
intensity after the extreme collimation (0.8 x10~° 
rad) is still rather strong: ~3x10° equivalent 
quanta/ minute. 

On the other hand we have determined the num- 
ber of pairs N(R, 6) as a function of k, for several 
fixed values of 6. 

Some of the results obtained are shown in this 
Letter. In Figs. 1(a) and 1(b) we give the ex- 
perimental ratio R,,: 


N(@,k) 01, (eo) 
N(6,k,) o,(k) , 





R,,(9)= (1) 


where o,(k)dk is the cross section for symmetri- 
cal pair production in aluminum at the photon en- 
ergy between k andk+dk. The solid line in both 
figures represents the value of the theoretical 
ratio Rp: 





In, @ quantity proportional to the bremsstrahlung 
intensity in the field of the nuclei of the single 
crystal, has been calculated by Uberall. 

Ig is a quantity proportional to the brems- 
strahlung intensity in the electron field of a non- 
crystalline target.° We have? 


I (8,x)= [1+ (1 - Ivo + v, (0/0) 


— A c 
+ B44 (0,8) - 80-aIfy,°C) 


+ ¥g(0/0)+ & vy"(@,0)], (a) 
h=1 


where x =k/E, E =1 Gev electron energy, and 
5 = (mc?/2E)[x/(1 -x)]=minimum momentum 
transferred to the nucleus in units of mc. 

The numerical values involved in formula (4) 
are the same as those used in reference 1, with 
the only difference that now we take the lattice 
spacing relative to the axis [111]. We computed 
the series 2 ~ 11", Lp - 12" only for the case 
6=0. 

In Fig. 1(a) we have k = 240 Mev and k,=910 
Mev. The experimental data show two symmetri- 
cal peaks at the left and right side of 0=0 anda 
central minimum in good qualitative agreement 
with Uberall’s calculations obtained in the Born 
approximation. Formula (4) results from an in- 
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FIG. 1. Intensity of the bremsstrahlung produced in a Si single crystal (T= 293°K) versus 0 (the angle between 
the 1~Gev electron beam and the crystal axis [111]). The solid curves of the figures represent Rip(9) given by for- 
mula (2), while the experimental points represent Rex() given by formula (1). The dashed curves are simply 
drawn for visualizing the behavior of the experimental points. The statistical error (~ 2%) is indicated for some 


points. (a) k=240 Mev; k,)=910 Mev. (b) k=80 Mev; k)=865 Mev. 
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effect among nuclei shows a strong dependence 
on the angle of the emitted photon.* This enables 
us to explain the fact that the maximum value of 
the measured effect is larger than the theoreti- 
cal one; this might be due to the enrichment 
suffered by the y-ray beam in photons emitted at 
small angles, owing to the sharp collimation em- 
ployed. 

In Fig. 1(b) we have k =80 Mev and k,=865 Mev. 
In the measurements shown in this figure the 
localization of the central minimum is very cri- 
tical. The spacing between the two maxima is 
about 2 mrad and the curve is very sharp, while 
the goniometric device has an angular repro- 
ducibility of + 0.5 mrad. We then proceed in the 
following manner. For small angles we have 


- 2 2) V2 
0=(8,7+0 2)", 


where 6), 9, are the angles of rotation of the 
crystal about a horizontal and a vertical axis, 
respectively. 

After a preliminary approximate alignment of 
the crystal axis with the electron beam, we ro- 
tate the crystal about a horizontal axis until we 
find a relative minimum, for which we put 6, =0. 
We then rotate it about the vertical axis by an 
angle 58,,=8,,*=1 mrad, to find the absolute 
minimum for which we put 0, = 0 (remember that 
the effect has an axial symmetry). 

In Fig. 1(b) we plot the values R,, obtained by 
rotation about the horizontal axis (crosses) ver- 
sus the angle 


=(@ 249 *2)12 
() (6, +8) , 


and the ones obtained by rotation about the ver- 
tical axis (circles) versus @=6,. 

Since the error 46/6 is very small for \0,/8y"| 
<<1, we have good angular resoiution for the 
crosses near the maxima. 

The analysis of the data plotted in Figs. 1(a) 
and 1(b) shows that the central minimum does 
exist, even if less marked than the theoretical 
one, especially as far as Fig. 1(b) is concerned. 
At this date we do not know if this difference be- 
tween theory and experiment really exists or if 
it is due to insufficient angular resolution of the 
experimental arrangement. A better goniometric 
device is under construction for a more precise 
analysis of the effect. 

At these energies it seems that the correction 
to the Born approximation near the minimum is 
smaller than the one predicted by Schiff,* accord- 
ing to which the central minimum would be com- 
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FIG. 2. Intensity of the bremsstrahlung produced in 
a Si single crystal (7 = 293°K) versus x =k/E (the frac- 
tional energy of the photon with respect to the electron 
energy E=1 Gev). The solid curve represents /(x, 6) 
given by formula (3) for 6=6 mrad. The dash-dot line 
represents the same quantity for a noncrystalline tar- 
get [H. Bethe and W. Heitler, Proc. Roy. Soc. (London) 
A146, 83 (1934); Wheeler and Lamb, reference 4]. The 
experimental points represent the quantity (5). The 
circles are relative to 8=6+0.5 mrad, while the 
crosses are relative to 6=1+0.5 mrad. 


pletely washed out for k=100 Mev. 

In Fig. 2 we plot the quantity J(x,6) given by 
formulas (3) and (4), as a function of x, for 6=6 
mrad (solid curve) and the bremsstrahlung inten- 
sity from a silicon noncrystalline target for com- 
parison (dot-dash curve). We also plot the ex- 
perimental quantity 


AN(k, 6)/o,(k), (5) 
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where A is a normalization factor, so chosen that 
the experimental value is coincident with the 
theoretical one at 900 Mev. The data for 6 mrad 
confirm that the collimation operates an angular 
selection of the photons. For 6~=1 mrad we do 
not draw the theoretical curve because at these 
angles the dependence on 6 is very strong, and 





for the experimental data we have 9=1+0.5 mrad. 


The preceding experiments” ® have shown a de- 
pendence on 6 of the bremsstrahlung intensity 
without showing any central minimum. We think 
that this is due to insufficient angular resolution 
with respect to the low-energy detected photons. 

We thank Dr. G. Barbiellini for the collabora- 
tion given in this work. 
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MUON ABSORPTION IN LIQUID HYDROGEN* 


Steven Weinberg 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received May 2, 1960) 


The long-awaited measurement of the rate of 
» absorption on protons could conceivably use 
hydrogen in any form, but there are technical 
advantages in working with liquids, as in bubble 
chambers. However, many physicists have ex- 
pressed pessimism about the utility of such an 
experiment, because in liquid H, muons form 
)-u -p molecular ions in addition to p- atoms. 
We shall show here how it is nevertheless possi- 
ble to interpret the observable rate of muon 
absorption in liquid hydrogen in terms of the 
tasic pnuv interaction. The main problem is to 
calculate the rate for the process 


(p-U-p)—-n+v+p (1) 


in terms of the muon absorption in atoms; be- 
cause of the extreme spin dependence of the 

V-A interaction,’ this requires knowledge of the 
orientation of proton and muon spins as well as 
the overlap of their wave functions at the instant 
of absorption. Almost all of our remarks will 
apply to solid as well as to liquid hydrogen, and 
we will confine ourselves to the case of isotopic- 
ally pure (H’),. 

The only two bound orbital states of the p-yu -p 
system” are a para lsog ground state with rota- 
tional angular momentum L =0 and binding en- 
ergy 2771 ev, and a 1sog ortho state with L =1 
which is AE =148 ev higher. The formation of 
these states by electron ejection in a collision of 













a p-wp atom with an H, molecule is respectively 
an E0 or E1 process, and hence the ratio of the 
formation rates is of order*® 


Para/Ortho ~ kv” ~2x10%, (2) 


where k, is the wave number of the ejected elec- 
tron. Actually, detailed calculations? show this 
ratio to be less than 3x10, so virtually all 
p-wu-p’s are formed in the ortho state. 

It seems furthermore that the p - u -p’s formed 
in the ortho state stay there during the few micro- 
seconds of muon lifetime. If we define the ortho— 
para electric dipole transition matrix element 
to be 


\(P | fpR10) | =ea,D, (3) 


then the rate for de-excitation accompanied by 
ejection of an electron from a hydrogen molecule 
in a collision is 


w (O-P)=162n m a *e*D?/3k 
e eeu e 


=2.7x10'°D? sec". (4) 


(Here n ¢ is the number density of electrons in a 
bubble chamber, n, =3.5x10" cm~*, and k, 

= (2m, (AE - Ey)”, where AE = 148 ev, and Ey 
= 15.6 ev is the electron separation energy in 
H,.) It is also possible that the p-u -p forms 
“ordinary” atoms or molecules with electrons 
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and protons from the liquid. The rate wa for 
Auger ortho ~ para conversion in such a system 
is again given by (4), with ~, replaced by the 
density x,’ of bound electrons at the p - u -p. 
If, for example, we take n,’= 1/na,*, then 


w ,(O~P) =61w (OP) =1.6x10"*D? sec™. (5) 


Of course, both w, and w, are much larger than 
the radiative conversion rate, 


w (O -P)= 4(AB)°a | %e*D*/S 


=8.1x10"D? sec™. (6) 


According to (4), (5), and (6), we would need 
D*2107 for any significant ortho — para conver - 
sion to occur during the muon lifetime, even if 
electrons were bound to all p-yu-p’s. It is un- 
likely that the £1 transition moment is this 
large. The ortho and para states have total pro- 
ton spins S, =1 and S,=0, and the only effects 
that can give rise to an E1 matrix element be- 
tween such states are impurities in the states,°® 
or the electric moments generated by moving 
magnetic dipoles.® The first effect would give 
at most a D of the order of the amplitude of the 
lsou (L =0, Sp =1) admixture in the para state, 
which can be estimated roughly as a typical 
magnetic interaction energy, a*m yp ~0.3 ev, 
divided by the minimum energy difference be- 
tween lsog and lsou states whose wave functions 
overlap appreciably, roughly 3 kev, giving an 
amplitude of about 10™*. Of course D is much 
smaller than this, since the lsou para impurity 
will not have perfect overlap with the lsog ortho 
wave function. The second effect mentioned 
would contribute to D an amount of order 
2.7T94E /m,~4x107". If D were this small, E1 
conversion would take as long as would be ex- 
pected for conversion by M2 transitions. 

De-excitation by rearrangement collisions of 
the p-yu -p with H, molecules is not retarded by 
the necessity to flip spins, but is hopelessly 
slow because of the high Coulomb barrier. The 
cross section for ortho — para conversion by this 
process contains the familiar factor eG where 
G =2n/1378 and 8c is the initial p -u -p velocity. 
In order that G<100, the p-yu -p would have to 
have kinetic energy greater than 200 ev, and this 
is impossible as it is formed with 0.02 ev, and 
kT <0.003 ev. 

The fact that process (1) occurs from the ortho 
state would not complicate matters badly if it 
were correct (as seems to have been previously 
assumed) that ortho p - up - p’s had definite total 
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spin 8, +8») =S(S+1). It has been shown’ that 
if w(S) is the rate for (1) in such states, and Wp 
is the rate for muon absorption in p- py atoms 
with total spin F, then 


w(1/2) = 2y (wp +40), (7 
w(3/2) =2y 5w,, (8) 


which may be compared with the smaller absorp- 
tion rate in the para state, 


w(para) = 2y p(4%> + q,). (9) 


Here yo, yp are the ratios of the muon density 
at one of the protons in an ortho or para p-u-p 
molecular ion (averaged over the separation be- 
tween protons) to the muon density at the proton 
ina p-yp atom. Using the CJR wave functions, 
we have found that® 


2y .=1.165, 2y_ =1.308. (10) 


O P 


These values are subject to corrections® of order 
my /mp due to admixture of higher orbitals than 
lsog. 

Of course, S is a good quantum number only 
if the forces coupling S» and Su to L are ignored. 
In this lowest approximation, the ortho state 
forms a degenerate quintuplet of states |J, S) 
which we can take with definite values of S and 
(Sy +Sp + L)* =J(J+1); the states are: 14, 3), 
1$,%), 18,8), 12,, 1%,%. The degeneracy 
among them is actually lifted by magnetic and 
relativistic effects which split the quintuplet into 
five states |y,) separated by energies of order 
a*m,, =0.3 ev, with total angular momentum 
J, =J,=43, Js=J,=%, J,;=§. [Each of these five 
states still retains a (2J,, +1)-fold degeneracy, 
which plays no role here. The field needed for 
a Paschen-Back effect is of order 10° gauss. | 

It is easy to show that the S= 4 and S = 3 com- 
ponents of the true energy eigenstates 1D») don’t 
interfere in reaction (1) if no attempt is made to 
measure the momentum or polarization of the 
left-over proton. Furthermore there is no inter- 
ference among the ld») states, since their en- 
ergy separation is much larger than their natural 
width of 3x107*° ev. Hence the rate for (1) is 
given as the sum 


= - 11 
up Ew(1/2) +(1 - €)w(3/2), (11) 
where 


b= Pte (12) 
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- = 2 = 
=I S=aly | ’ (§, =0) (13) 
and P, is the probability of the ortho p-yu -p 
occupying lp, at the instant of absorption. 
It is easy to calculate the relative probabilities 
Ny for formation of the five ortho states. Since 
the p -u -p is formed from a proton and an 


F=0 p-u atom,” it has S=4 initially, and hence 
Se. - 
P = 42d +1)E (14) 


Furthermore, it seems reasonable to suppose 
that no transitions among the |y,) occur during 
the muon lifetime. The Coulomb barrier pre- 
vents rearrangement collisions’; M1 radiative 
transitions would take a day; and probably M1 
conversion in collisions with H, molecules is 
very slow because the energy available is in- 
sufficient for electron ejection. Assuming then 
that P, =P,°, we obtain 


E=D 42d, +0E,?, (15) 


and our problem reduces to that of computing 
the £,. However, there is one important in- 
equality that can be derived without any calcula- 
tion. Since £, += +,=1, it follows directly 
from (15) that 


5<é<1, (16) 


whereas in general we would only know that 
0<&<1. If we accept that w,~50w,,” then be- 
cause we cannot have é near zero it is safe to 
approximate 


x= /W9 = (27 5) 48, (17) 


Ww 
pup 
and furthermore the neutrons emitted should 
have almost 100% longitudinal polarization. 
__ The Hamiltonian H’ that couples Sp; Ss we and 
L takes the general form 
H’ a _ » ~_ 
E 15» L+E 2 s°p p 
+E [XS -L)S.-L)+3(8 -L 
4l26 MS, ) a(S, 
+E_[(S, -L)?+3(8. -L) - 4]. (18) 


rs ~-L+E $s 8 
yb 


MB h)- 56, °8)] 


One term in E, represents the proton spin-orbit 
force; E, and another term in E, come from the 
magnetic interaction between the proton momen- 
ta and the muon and proton spins; E, and E, 

come from the interaction between proton and 
muon magnetic moments; and E, arises from the 
interaction between the proton magnetic moments. 
We have derived formulas (too lengthy for this 


note) which give the é, as simple algebraic func- 
tions of the E j coefficients, and which express 
the E; in terms of averages over the molecular 
ion orbital wave function. All the E; seem rough - 
ly of the same magnitude, containing factors 
e*/mymy or e?/my’s e.g., 

Ey = Be*(2.79)y ,/3m ma ° =0.292 ev. 
If E, -E, and E, -2E, turn out to be small com- 
pared to E,, then — will be near unity. We have 
not performed the necessary integrations to 
obtain the E; and hence ~ numerically, as it 
seems most important at this stage to recognize 
that the task is straightforward. 

Since some muon absorptions will occur with 
muons still in F =0 p-u atoms, the observed 
total rate of all absorptions in liquid hydrogen 
iS Waps = Wo(1+xy)/(1+y), where y is the product 
of the muon mean life times the rate per p-wu 
atom of p--p formation. According to CJR, 

y =14.4 at a proton number density n, =3.5 x10”? 
cm’, Since x [as given by (17), (10), and (16)] 
must be between 0.44 and 0.88, the value of w, 
determined by a measurement of Wapg (and cal- 
culation of x) will be quite insensitive to the 
precise value of y, providing y isn’t too small. 

It is possible to test many of our assumptions 
by measuring w,p, Over a range of hydrogen 
densities. If CJR are correct, varying Np from 
1.8 to 4.5 x10” cm™ (possible in a counter ex- 
periment) would change y from 7.4 to 18.9, 
giving according to (17) a decrease in wy, of 
between 1% and 8%. If an earlier calculation,’ 
giving a y five times smaller than that of CJR, 
were correct, then the decrease in Waps would 
be between 2.5% and 16.5%. If our estimates 
are wrong, and ortho — para conversion becomes 
important at the higher densities, x and wap, 
can drop by almost a factor of 3, while if ortho — 
ortho transitions became important wap, could 
rise or fall. 

The author has benefited from conversations 
on these matters with many physicists, and 
would particularly like to thank R. H. Dalitz, 

R. Hildebrand, D. L. Judd, L. M. Lederman, 

R. J. Myers, W. Nierenberg, H. Primakoff, 

M. A. Ruderman, R. Sard, J. Steinberger, V. L. 
Telegdi, and S. B. Treiman. 
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P-WAVE RESONANCE IN PION-PION SCATTERING* 


Benjamin W. Leet and Michael T. Vaughn 
Department of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received April 25, 1960) 


The possible p-wave resonance in pion-pion 
scattering has been a subject of much discussion 
lately.‘“* The work of Chew and Mandelstam 
based on the double dispersion representation® 
for the transition amplitude suggests that there 
is more than one solution to the set of integral 
equations for the pion-pion partial wave scatter- 
ing amplitudes.*** In the present note, we wish 
to point out the possible source of such a multi- 
plicity of solutions from the viewpoint of conven- 
tional Lagrangian field theory, and to examine 
the significance of the fact that the p-wave phase 
shift remains small throughout the physically 
acceptable range of the pion-pion coupling con- 
stant in the “s-wave dominant solution” of the 
Chew-Mandelstam equations,*** although a solu- 
tion containing a p-wave resonance may be ob- 
tainable from a nonadiabatic approach, the start- 
ing point of which we suggest here. 

It is crucial for the following discussion to 
realize that the form of the double dispersion 
representation® is uniquely determined by the 
masses of the participating particles and the se- 
lection rules, aside from the subtle question of 
possible subtractions necessary to give the re- 
presentation a valid meaning. If one assumes 
the existence of a vector boson of isotopic spin 
one, which interacts with the isotopic vector 
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part of the pion current in the Lagrangian,® and 
insists that the renormalized mass (which we 
will define later) of the postulated particle be 
larger than twice the pion mass, then the par- 
ticle will become unstable and one has not changed 
the selection rules of the theory from those of the 
conventional pseudoscalar meson theory. Sucha 
theory, on the other hand, will necessarily pre- 
dict a p-wave resonance in pion-pion scattering. 
The situation here is similar to the well-known 
Castillejo-Dalitz-Dyson ambiguity’ in the static 
meson theory: the double dispersion relation 
does not imply one particular Lagrangian, but a 
class of Lagrangians which are selection-rule 
equivalent.® 

It follows then that there exists a solution to 
the Chew-Mandelstam equations corresponding 
to a Lagrangian theory in which there exists, in 
addition to the usual couplings, a coupling term 
of the form 





1 i j ij k 
= .a(p 8 -8 B (1) 
where B, * is the field operator for the vector 
boson of isotopic spin one, & being the isotopic 
spin index, and ¢ is the pion field operator. 
The p-wave resonance in pion-pion scattering 
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is then due to the resonant structure of the vec- 
tor boson propagator, the relevant part of which 
is 


pe) = -D(R*) g*” (2) 


(the terms proportional to k¥kY do not contribute 
to pion-pion scattering in the present approxima- 
tion), where 

D(k?) = [k? - m? - Z(k?)] 3, (3) 


where m>2y is the renormalized vector boson 
mass, and D(k*) is the vector boson self-energy 
operator; in the “resonance approximation,” in 
which we consider only the sum of the iterated 
bubble diagrams (see Fig. 1), D(%?) is given 
simply by 


5(k*) 
; LV, 
4i kk 

q 4 “~ ) 


v 
* 3(2n)4 i [q?- uw? | (@+k)?- v7] ” Re 





(4) 


The imaginary part of D(%”) can be calculated 
directly; it is given by 


Im 5(s +ie)=- 2 S-40™ 


J Se — 06-4u7). (5) 


The real part can be calculated, for example by 
noting that D(s) satisfies the dispersion relation, 


L(s) = = (m?) + (s-m?) D’(m?) 


lim @-m?) (ImZ(’) 
n-m* s’-m?*-ie 


x r A r 7 A hae, (6) 
S'-S-ie€ $'-n-ie 


where we have made two subractions at s =m?, 
reflecting the fact that the integral in Eq. (4) is 
quadratically divergent. Note that while £(m”), 
i'(m?) are now complex, their imaginary parts 
explicitly appearing on the right-hand side of (6) 
are cancelled by the contributions from the ie’s 
inthe denominators (s’-m?-ie)~' and (s’-n-ie)~; 
the real parts of these constants are absorbed 
inthe definition of the renormalized mass and 
coupling constant. Carrying out the indicated 
integrations then leads to an expression for the 
renormalized self-energy operator: 


Z(s) =J(s) - J(m?)+ilmZ(s), (7) 








_ f? S(s-4u?)™ | | s+ (s-4u?)” 
Js) = fea fe ee] . | (8) 
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FIG. 1. Diagrams included in resonance approxima- 
tion to boson propagator. 


and the constant é is given by 
_ m*-4y? 


m2 


m?+2u2\ (m?-4y2)¥? = m+(m?-4 2)" 
*1 on m m-(m?-42)|° 





é 





(9) 


The J=1, J=1 pion-pion scattering phase shift 
6 is then given in this approximation by 


V2 . 
4 Ps one 
1 


7 v Kv )-vK(v)-i r[/(v+1)]*’ 








(10) 


where v= (s-4y?)/4y?, v,=(m?-4u?)/4y?, D=f?/ 
487, and 
(v+1)*+ py” 
K(v) =1- EJ" In 1) pe . 


nm WN\v+l1 


We note that this expression for the scattering 
amplitude has a pole for v=- y, on the negative 
real axis; indeed, Eq. (10) was first written 
down, in a slightly altered form,® by Frazer and 
Fulco,’ who derived it from the partial wave dis- 
persion relations® with the left-hand branch cut 
replaced by a phenomenological pole. From the 
present point of view, this pole merely repre- 
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sents the failure of the resonance approximation 
for large |s|; in the resonance region, the effect 
of the pole is negligible if v, is sufficiently large 
and f is not too large. The distinction between 
the two viewpoints may be seen by noting that if 
the pole is interpreted as a phenomenological 
representation of the left-hand branch cut, then 
v,~m?, where m is some average mass of the 
intermediate states which contribute to the left- 
hand branch cut. Since v,~100 - 1000 for values 
of the width [ which fit the data on nucleon elec- 
tromagnetic structure,’ one is led to suspect that 
virtual baryon-antibaryon pairs play a prominent 
role in producing the p-wave resonance,*° un- 
less an unstable vector boson such as Sakurai 
and others have suggested® is present in the ori- 
ginal field-theoretical Lagrangian. 

If Eq. (10), with the pole removed, is inserted 
into the iterative procedure of Chew et al.,* then 
the resonance should persist, and, in addition, 
nonresonant terms will appear in other angular 
momentum states; in particular, one should be 
able to calculate the s-wave scattering lengths 
in terms of the resonance parameters. It also 
seems clear that a purely dispersion theoretic 
approach will require at least two parameters 
to obtain a p-wave resonance, since the one- 
parameter theory of Chew et al., which corre- 
sponds to a \¢* pion-pion interaction,'! seems 
incapable of producing a resonance. Whether 
these two parameters are to be related to the 
pion-pion and pion-nucleon coupling constants, 





or to the mass and coupling constant of the con- 
jectured vector boson, cannot be decided on the 
basis of present knowledge. 





*Supported in part by the U. S. Atomic Energy Con- 
mission. 
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"L. Castillejo, R. Dalitz, and F. Dyson, Phys. Rev. 
101, 453 (1956). 

®This point is discussed for the static theory by R. 
Norton and A. Klein, Phys. Rev. 109, 991 (1958). 

More precisely, I’ = nI(7+éT) = vy’ = I’ vy K (vy), 
where vy’, I’ are the quantities designated as v,;, Ty 
by Frazer and Fulco.'*? 

1K, Igi and K. Kawarabayashi, Progr. Theoret. 
Phys. (Kyoto) 20, 576 (1958), have noted that the con- 
tribution to the P-wave pion-pion scattering amplitude 
from a single nucleon closed loop corresponds to an 
attractive potential. 

11§, Okubo, Phys. Rev. 118, 357 (1960). 





NEW THEORETICAL VALUE FOR THE LAMB SHIFT* 


Arthur J. Layzer 
Columbia University, New York, New York 
(Received May 9, 1960) 


Although the good agreement between theory 
and experiment for the Lamb shift is one of the 
notable successes of quantum electrodynamics, 
the most recent tabulations’ still have shown 
small but significant discrepancies for H, D, 
and He* which have made a further increase in 
the accuracy of the theoretical value desirable. 
Since the listed discrepancy for He” is roughly 
60 times that of H or D, it was natural to try to 
reduce this disagreement by calculating the 
second order radiative correction of order 
a(aZ)*, which is the next order in(aZ), the 
Coulomb interaction parameter, beyond the pre- 


580 


viously calculated’»* order of a(aZ)°. 

A closer mathematical analysis*»® shows that 
there exist the leading orders of a(aZ)*In*(aZ) 
and a(aZ)*ln(aZ). We have completed the cal- 
culation of these two orders. The analytic result 
is 

AE(2S - 2P,,) 


= -Lw{ } In?w + Inw(4 In2 + 1+ 7/48)], (1) 


where AE(2S -2P,,.) is the difference of the shifts 
of the 2S and 2P,, levels due to the two new orders; 
w=(aZ)*, and L is Z* times the “Lamb constant”. 
L=Z*a°/(3n) ry =Z*(135.6) Mc/sec. The aw*ln’w 
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coefficient vanishes for non-S states and is pro- 
portional to 1/n* for the nS state. The value of 
the aw°ln*w coefficient has been recently con- 
firmed by an independent calculation of Fried 
and Yennie.°® 

The corresponding numerical results in Mc/sec 
are -0.25 for H or D and -9.5 for He*. These 
additions reduce markedly the disagreement be- 
tween theory and experiment and, when added to 
the theoretical values listed by Petermann,’ lead 
to the new theoretical values listed in Table I, 
which are compared with experimental values 
for H, D, and He*. Theoretical errors are copied 
from Petermann’s article. 

The theoretical foundation of the calculation of 
the new orders has been briefly described in 
reference 4. The method used is the free- 
propagator expansion, the algebraic expansion 
of the bound electron propagator or Green’s 
function in “powers” of the Coulomb potential. 
This expansion was first employed in the Lamb 
shift problem by Fried and Yennie*® who success- 
fully extracted the two lowest orders of Linw 
and L with the aid of a special photon gauge. In 
our method, this gauge is not employed and the 
photon propagator is of the usual Feynman form 
proportional to 1/k’?. 

Details of the calculation of the two new orders 
will be published later along with more general 
results of a mathematical nature concerning some 
properties of the free-propagator expansion con- 
nected with an expansion of the self-energy in 
orders of (@Z). The remainder of this report 
is devoted to a further presentation of the re- 
sults of the calculation for arbitrary bound states 
and in particular for the 1S, 2S, 2P,,., and 2P,, 
states. 

Let us introduce the notation 


¥7(0) = | fd8p{ piu) |? =8n3 | (x =Olu) I?, (2) 
¥?(0) (mn?) =8/(n*) for nS state 
=0 for L #0, (3) 


where u is the nonrelativistic Schrodinger wave 


function in “Z-atomic” units (8 =c=aZm=1). 

These units will be understood in the following. 
The shift AE, due to both the aw*ln’w and 

aw*lnw terms can be conveniently expressed as 


AE =AE, +AE,, (4): 


where AE, is of order aw*lnw and proportional 
to ¥7(0), and AE, is the sum of a number of 
matrix elements which can be easily evaluated 
to the desired order in w for the particular bound 
state under consideration. AE, contains the en- 
tire aw*ln*w coefficient, proportional to 77(0), 
and the remainder of the aw*lnw coefficient, 
which is not proportional to (0). 

AE, is given by 

y?(0) 


AE, = -L—,—winw(4 In2 - % - 4=-0.26). (5) 


The last term in the parenthesis, -1/10, is due 
to vacuum polarization and can be easily derived 
as the effect of the Dirac modification of the 
large-component wave function on the expecta- 
tion value of the Uehling potential. 

The remainder of AE is given by (V = -1/r) 


az, --Lwf “2 
x[-4(u’!...¢?V lu) - (0, bP) 
- (11/10)(p,...4°V...D ) + (8/5){p?...q?V) 
+2(p*...[b, Vp) -6(L>,,V]...[b,,V])]- (6) 


Here uw’ is the first “Dirac correction” to the 
expansion in powers of w of the exact large- 
component wave function y: y=u+wu’. The 
angular brackets ( ) denote the expectation value 
with respect to the nonrelativistic Schrodinger 
wave functions u. The three dots denote a con- 
vergence factor of (1+wp?)"'. The lower limit 
of integration is arbitrary but finite and fixed 
(w- independent).” 

The evaluation of AE, for the 1S, 2S, 2P,,, 
and 2P,, states yields, with the omission of the 


Table I. Theoretical and experimental values of the Lamb shift in Mc/sec for H, D, and He*. 











H p* He* > 
Theoretical 1057.70 + 0.15 1058.96 + 0.16 14046.3 + 3.0 
Experimental 1057.77 + 0.10 1059.00 + 0.10 14040.244.5 


_— 





“8. Triebwasser, E. Dayhoff, and W. Lamb, Phys. Rev. 89, 98 (1953). 
E. Lipworth and R. Novick, Phys. Rev. 108, 1434 (1957). 
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previously given aw*ln*w term, 
AE,(1S) = -Lwlnw(24 In2 + 86/5 =33.8), (7) 
AE, (2S) = -Lwinw(369 /80 = 4.61), (8) 
AE ,(2P,/2) = -Lwinw(103 /240 = 0.43), (9) 
ME, (2P 12 - 2P 2) = -Lwinw(-3/16 =-0.19). (10) 


The contribution to the fine structure separa- 
tion 2P,.-2P,,, is due, of course, to the spin- 
orbit matrix element of (6). 

We conclude this report with a brief specula- 
tion about the magnitude of the uncomputed term 
of order wL, which is technically of higher order 
than the orders calculated here through the ab- 
sence of factors of lnw. Since the coefficient of 
the lnw term in the 2S -2P,,. shift given in (1) is 
about 5 times as large as the In*w coefficient 
and (-Inw) is only about 10 (84) for H (He), the 
contribution of the Inw term is roughly 1/2 (3/5) 
that of the In’w term for H (He*) with opposite 
sign. A comparable ratio of the “constant” to 
the Inw contribution would mean that the term of 
order wL is still important. 

There are theoretical reasons, however, for 
believing that this is too dark a picture and that 
the “constant” term is only about 10% (20%) of 
the shift calculated here for H (He*) due to the 
logarithmic terms. Let us note first that the ln? 
and ln contributions will change if one picks 
w’=,w rather than w as the argument of the 
logarithmic factors, where A is an arbitrary 
w-independent number, a substitution which cor- 
responds to the choice of w’ rather than w as the 
“natural” expansion parameter.* Thus, Aln*w 
+ Blnw + C =Aln?w’ + B(A)lnw’+C(a), where B(a) 
= B-2Alna and C(A) =C +Aln’a - Bina. 

If one assumes that the uncalculated “constant” 
term, C(A), is negative for A=1, which is in 
accordance with the fact that all five previous 
orders in w alternate in sign for A =1 (and w<1), 
then it is rigorously true that the “constant” term 
is of minimum absolute value (and <0) for some 
value of A greater than unity but less than A* 
where A* is such that the coefficient, B(A), of 


the Inw’ term vanishes. In our case, A* =e?** =13. 


With the choice A =A* one finds that the contribu- 
tion in Mc/sec to the Lamb shift due to the 


logarithmic terms is changed from -0.25 to -0.% 
and from -9.5 to -11.7 for H and He*, respec- 
tively.*® It is satisfying that the change is small 
and in a direction towards improving the over-al| 
agreement with experiment. 

I would like to express my appreciation to 
Professor N. M. Kroll for suggesting the problem 
of calculating the new orders and for much help. 
ful advice and encouragement. I wish to thank 
Dr. D. R. Yennie and Dr. H. Fried for several 
interesting discussions. 





This work was supported in part by the U. S. Atomic 
Energy Commission. 

1A, Petermann, Fortschr. Physik 6, 505 (1958). We 
refer to this article also for a discussion and classifi- 
cation of the numerous theoretical effects. 

2M. Baranger, H. A. Bethe, and R. P. Feynman, 
Phys. Rev. 92, 482 (1953). 

5R. Karplus, A. Klein, and J. Schwinger, Phys. Rev, 
86, 876 (1952). 

A. Layzer, Bull. Am. Phys. Soc. 4, 280 (1959). 

5H. M. Fried and D. R. Yennie, Phys. Rev. 112, 
1391 (1958). 

6H. M. Fried and D. R. Yennie, following Letter 
[Phys. Rev. Letters 4, 583 (1960)]. 

"The spin-orbit and (p;j...q?V...p;) matrix elements 
vanish for S states. The Dirac and (p?...q?V) matrix 
elements vanish for non-S states. For non-S states the 
convergence factor may be eliminated and { Wow /w re- 
placed by Inw. For S states, the nonvanishing Schrod- 
inger matrix elements diverge like w~”? in the limit 
w-— 0 and therefore give rise to a lower order aw” 
contribution to AE that should be discarded. [The 
w~”? part of the matrix elements could be eliminated 
by replacing w by -w+ie in the convergence factors 
(1+wp*)~! and taking the real part of the matrix ele- 
ments.] The aw*ln*w contribution is due to matrix 
elements with a Inw divergence for S states. This 
divergence is proportional to ¥7(0)._ The corresponding 
matrix elements are the Dirac matrix element and the 
two Schrodinger matrix elements involving [);, V). 

8A similar situation arises for the lowest order con- 
tributions of order aw*Inw (~Linw) and aw? (~L). 
There, as is well-known, it is natural to introduce an 
“average excitation energy’ for the 2S state of magni- 
tude (17. 7)/2=8.8 and the term of order aw’? can be 
made to vanish for a choice of A nearly equal to this 
excitation energy. 

*Virtually the same results would have been obtained 
for the choice A= 8.8 corresponding to the 2S average 
excitation energy. 
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Following previously described techniques,’ 
the atomic energy level displacements of order 


a result in agreement with that of the preceding 
Letter.? The procedure utilized represents a 
direct extension of the methods introduced in 


having zero expectation value for the particular 
atomic state in question, is added to the series 


state electron propagator, >) - I,- This addi- 
tion corresponds to a change of gauge of the vir- 
tual photon defining the one-photon Lamb shift, 
and serves to remove spurious lower order con- 
tributions from each term of the combined se- 


calculational scheme for the lowest order Lamb 
shift and a particularly convenient method for 
the estimation of the magnitude of the Bethe 
logarithm. 

It may also be of some interest to remark that 
the first relativistic corrections to the lowest 
order Lamb shift [of order a(Za)*mc?*], and to 
hfs [of order a(Z aE ys), may both be obtained in 
a quite straightforward manner by the use of 
these techniques. In fact, the necessary matrix 
elements for the hfs calculation can be read off 
immediately from those of the similar Lamb 
shift effect. It is easily seen, for either case, 
that contributions of the desired order arise 
only from the terms M,+M,+M,. A simple 
parametric decomposition uniquely isolates that 
portion of these matrix elements of the proper 
order (in Za); and the evaluation of elementary 
integrals then yields the correct, previously ob- 
tained® results. 

Corrections of order a(Za)® In*(Za)mc? to the 
Lamb shift* arise from the matrix elements of 
M,+M,+M,+M,, where the contributions of M, 
are all, to this order, gauge-variant, and serve 
to cancel similar gauge-variant terms remaining 
in the combination M,+M,+M,. As a convenient 
check on the manipulations, all matrix elements 
may be grouped into gauge-invariant and gauge- 








a(Za)® In?(Za)mc? have been calculated and yield 


reference 1, wherein a counter-series, Linz dy 


representing the iterative expansion of the bound- 


quence, M,, =I,+J,, thereby generating a simple 





HIGHER ORDER TERMS IN THE LAMB SHIFT CALCULATION 
H. M. Fried 
University of California, Los Angeles, California 
and 


D. R. Yennie 
University of Minnesota, Minneapolis, Minnesota 
(Received May 16, 1960) 


variant combinations, and the latter explicitly 
shown to cancel. To this order in Za, it is 
necessary to include a relativistic correction to 
the ordinary S-state atomic wave functions bn”); 
this essentially corresponds to averaging the 
latter, together with that portion of the operator 
which leads to the Bethe logarithm in the lowest 
order calculation, over distances comparable to 
the electron’s Compton wavelength. To illustrate, 
consider the large component of the 2S-state 
(Dirac) wave function which, correct to the or- 
der a(Za)® In*(Zq) in the final result, may be 
written in the form 


(7) = | o,(0) ie (1 -yryr)° sa 


where s =[1 -(Za)*]”, and y=Zam/2. Approxi- 
mating the quantity (yr)S~1 by 1-4(Za)? In(y7), 
one obtains the desired correction 


4¢,(r) =-3(Za)? Inr)¢,(7). 


In conjunction with that part of M, which yields 
the Bethe logarithm (corresponding to the infra- 
red divergent part of the free-particle vertex 
operator), this wave function correction (wfc) 
will yield a relative factor of (Za)* In(Zq@), since 
the corresponding matrix elements will contri- 
bute only when r~m~! (momentum ~m, rather 
than Zam). 

However, the major contributions of this or- 
der arise from the structure of the operators 
themselves; in essence one has to deal with in- 
tegrals of the form® 


1,17 4 ‘dex oo)", (1) 
0 


where the index specifies the parent matrix 
element, viz.: 


D,@) =xm?+(1 -*)A, +xnry(1 - y)q? 
appears in the terms of M,. The parameter \ 
may have any positive value, including zero; 
and {dT stands for some combination of the di- 
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mensionless momenta, obtained from those ma- 
trix elements of the correct nominal order 
(without wfc), which would be logarithmically 
divergent without the convergence factors D, (a). 
Omitting that part of M, which leads to the wfc 
contribution, the gauge-invariant terms which 
appear in M,+M, may be written as 


2(),xp)?a + 2imaalg P* - +Q), 
while those of M, are 
2imaa(- 2 P* - 2 @), 


where P, Q denote the “gauge-invariant” mo- 
mentum combinations 

P= 2p"-p,?-p, = G,-p’)?+ @-b’)? - @,-D”. 
The integrations over momenta and the para- 
metric variables, as indicated in Eq. (1), have 
been suppressed; it is found that the Iln*(Za) con- 
tributions of such integrals are independent of 
l,m, and the denominator “y variables.” These 
contributions may conveniently be determined by 
noting that if one anticipates Za dependence of 
the form 


I(Z a) =A In*(Za) + Bln(Za)+C+D(Za)+..., 


the desired coefficient, A, may be obtained by 
calculating 


a 8 al(Z a) 
2 [ea weai(20) pl —, 


The resulting integrals are vastly simpler to 
evaluate than are the original J(Zq@); and in this 
way it is not difficult to obtain the final answer 
(nonzero for S states only): 


AE = a a(Z a)® ln*(Z a)me{1 - 4}. (2) 

The first term in the bracket of Eq. (2) repre- 
sents the wfc contribution; the second one arises 
from elements of the self-energy operator, as 
outlined above. The numerical effect of this re- 
sult on the theoretical estimate of the Lamb 
shift for H, D, and Het is presented in the pre- 
ceding Letter. 

We wish to thank Professor E. Wichmann, 
Professor N. M. Kroll, Professor H. Suura, and 
Mr. A. J. Layzer for interesting conversations 
relative to this work. 





‘H. M. Fried and D. R. Yennie, Phys. Rev. 112, 
1391 (1958). 

2a. Layzer, preceding Letter [Phys. Rev. Letters 4, 
580 (1960)]. - 

3R. Karplus, A. Klein, and J. Schwinger, Phys. 
Rev. 86, 288 (1952); M. Baranger, H. A. Bethe, and 
R. P. Feynman, Phys. Rev. 92, 482 (1953); N. M. 
Kroll and F. Pollock, Phys. Rev. 86, 876 (1952); R. 
Karplus, A. Klein, and J. Schwinger, Phys. Rev. 84, 
597 (1951). 

‘The probable existence of such corrections was 
pointed out several years ago by E. Wichmann (private 
communication). 

5We use the notation of reference 1. 
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NEW TEST FOR Al=3 IN K* DECAY. Steven 
Weinberg [Phys. Rev. Letters 4, 87 (1960) ]. 


It was incorrectly stated in Eq. (17) that the 
parameters aj and y; characterizing the energy 
distribution of the unlike and like pions in 7 or 
7’ decay are related, as a consequence of 
Eq. (14), by Y¥j= -a;/2. Dr. R. G. Glasser has 
kindly informed me that the relation should read 


-a. 
J 
ae ree ya (17)’ 
j 2+ a, @; /m K 
and that use of this correct result brings theory 
and experiment into remarkable accord. [Equa- 
tion (17)’ holds strictly only for j=7.] The ex- 
perimental value (reference 6) of y; is -2.2+0.3, 
which may be converted using (17)’ into an equi- 
valent a, =6.6+0.9, whereas the directly meas- 
ured value is a-=6.8+1.2. This close agreement 
indicates that Eq. (14) is an excellent approxi- 
mation, all departures from a phase space dis- 


tribution being accounted for by the linear factor 
1+fQy/mx, and hence that the 7’ spectrum may 
be predicted with confidence. From the above 
values of a, we obtain 8,=1.65 (agreeing with 
Dalitz’s value quoted in footnote 7 as 8, =1.6+ 0.5) 
and predict that 8, = -3.3 and hence that a_, = -6.3. 
A recent analysis of 72 7’ events’ gives a value 
a_,=-7.1+3, so that AI =3 seems to be working 
well. 


1S. Bjorklund, E. L. Koller, and S. Taylor, Phys. 
Rev. Letters 4, 424, 475(E) (1960). 


PROPOSAL FOR AN ELECTRON SPIN RESO- 
NANCE EXPERIMENT OF S-STATE IONS 
UNDER HIGH HYDROSTATIC PRESSURE. 
Hiroshi Watanabe [Phys. Rev. Letters 4, 410 
(1960) ]. 


The fraction (1/336) in the expression (2) 
should read (1/168). 
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ABSTRACTS 








In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


SOLUBLE THREE-DIMENSIONAL MODEL FOR 
TOWNSEND’S a. G. W. Stuart and E. Gerjuoy, 
John Jay Hopkins Laboratory for Pure and Ap- 
plied Science, General Atomic, San Diego, Cali- 
fornia (Received March 18, 1960). 


A model gas is considered in which all elec- 
tron-molecule cross sections are isotropic and 
depend inversely on the velocity v. Collisional 
energy loss is neglected. The Boltzmann equa- 
tion for the model is solved for the collision 
density, where the collision density is the num- 
ber of collisions that an individual electron 
makes between v and v+dv over its entire his- 
tory. The Townsend a is obtained from the col- 
lision density, and it is found that a/p is in- 
versely proportional to E/p. It is argued that 
this model furnishes an upper bound to the true 
a/p for all E/p; therefore it is concluded that 
this model demonstrates that at sufficiently high 
E/p the observed a/p for any real gas must de- 
crease with increasing E/p. The results also 
shed light on the way electron energy balance or 
lack of energy balance affects a/p and the drift 
velocity vp; it is shown that energy balance is 
not possible at arbitrarily large E/p. Numerical 
applications of these results to H, are discussed. 


STUDIES OF PLASMA HEATED IN A FAST- 
RISING AXIAL MAGNETIC FIELD (SCYLLA). 

K. Boyer, W. C. Elmore, E. M. Little, W. E. 
Quinn, and J. L. Tuck, Los Alamos Scientific 
Laboratory, University of California, Los 
Alamos, New Mexico (Received March 28, 1960). 


The Scylla plasma experiment, which employs 
a rapidly rising magnetic field in a cylindrical 
mirror geometry to produce and heat a deuterium 
plasma, is described. Experimental studies of 
the reproducible neutron emission from the hot 
plasma show that the neutrons are emitted (1) in 
a symmetrical, bell-shaped time distribution 
centered on the maximum of the magnetic field, 
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(2) from a limited region with a 2-cm axial length 
and a 1.5-cm diameter centered in the compres- 
sion coil, and (3) in the radial direction with a 
narrow spread of energies and no significant 
anisotropy. The time distribution of the neutron 
emission is shown to be in agreement with a 
thermonuclear yield curve calculated for an 
adiabatic compression by the observed magnetic 
field. The neutron yield has been studied as a 
function of deuterium pressure, capacitor -bank 
voltage, and nitrogen impurity. Observations of 
the space-time distribution of the visible light 
emission with a streak camera show that (1) a 
strong radial “ shock” occurs at the beginning of 
the second half-cycle, (2) very little light is 
emitted from the plasma “fireball” during the 
time of neutron emission, and (3) an intense 
luminous flux is produced during the later stages 
of the discharge. The energy absorbed in each 
half-cycle of the discharge by the gas is presented 
as calculated from the incremental damping of 
the driving magnetic field. Observations of hard 
x-ray emission (~ 200 kev) at times of maximum 
dB/dt for operating pressures in the 5 to 50 
micron range are contrasted with the character- 
istics of the neutron emission in regard to time 
distribution, pressure, impurities, and rf pre- 
excitation. Magnetic probe studies of the Scylla 
discharge are reported and evidence is given 
that the perturbing effects of the probe dominate 
the plasma temperature. 








CONTINUUM RADIATION IN THE X-RAY AND 
VISIBLE REGIONS FROM A MAGNETICALLY 
COMPRESSED PLASMA (SCYLLA). F. C. Jahoda, 
E. M. Little, W. E. Quinn, G. A. Sawyer, and 

T. F. Stratton, Los Alamos Scientific Laboratory, 
University of California, Los Alamos, New 
Mexico (Received March 28, 1960). 


The identification of a sharp low-wavelength 
cutoff in the spectrum of x rays emitted from 
deuterium discharges in Scylla has resulted in 
the assignment of an electron temperature of 
240+ 40 ev at the time of peak magnetic field 
compression. Simultaneous time-resolved abso- 
lute intensity determinations in the visible con- 
tinuum, when coupled with the temperature meas- 
urement, yield an upper limit electron number 
density of (5+ 1)x10'*/cm* at peak compression. 
The absolute value of dE /da in the soft x-ray 
region is two hundred times larger than brems- 
strahlung from a pure deuterium plasma at the 
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temperature and density quoted, and it is pos- 
tulated that the large experimental dE /dd is the 
result of recombination radiation from about 2% 
of oxygen contaminant from the discharge tube 
walls. 


VELOCITY SPECTRUM OF PROTONS AND TRI- 
TONS FROM THE d -d REACTION IN SCYLLA. 
D. E. Nagle, W. E. Quinn, F. L. Ribe, and W. B. 
Riesenfeld, Los Alamos Scientific Laboratory, 
University of California, Los Alamos, New 
Mexico (Received March 28, 1960). 


A diagnostic experiment has been carried out 
on the d -d reactions produced by fast magnetic 
compression of a deuterium plasma. A deter- 
mination of the velocity spectra of protons and 
of tritons from the d -d reaction was made by 
magnetic analysis and nuclear plate detection of 
the particles. The observed distributions are 
Gaussian, with widths which correspond to a 
deuteron temperature of 1.3 kev. Comparison 
of the mean proton and triton momenta indicates 
that no plasma drift in the (axial) direction of 
observation is present, nor any potential differ - 
ence between the source plasma and detector 
greater than a few volts. These results, coupled 
with previous ones on the neutron yield, duration, 
source extent, and lack of circumferential drift 
argue against any of the simple, physically plau- 
sible non-Maxwellian acceleration mechanisms 
for the d -d reactions so far proposed. 


DIFFUSIVE MOTIONS IN WATER AND COLD 
NEUTRON SCATTERING. K. S. Singwi* and 
Alf Sj6lander, Argonne National Laboratory, 
Lemont, Illinois (Received March 17, 1960). 


Using a model of liquid water in which a mole- 
tule, in its equilibrium position, performs an 
oscillatory motion for a mean time 7,, and then 
diffuses by continuous motion for a mean time 
t, and repeats this sort of motion, the differ - 
ential scattering cross section for cold neutrons 
tas been calculated. It is found that the shape 
ofthe “quasi-elastic” scattering is, in general, 
tot Lorentzian. The formula for the broadening 
of the quasi-elastic peak assumes a simple 
form in the two limiting cases: In case (i), 
'\T), it reduces to the formula derived on the 
simple diffusion theory; and in case (ii), T,«7,, 
the broadening is the same as in case (i) if 
y Dr,«1, and it approaches the asymptotic 


value 24/1, if x*DT,>1, where hx is the mo- 
mentum transferred to the system and D is the 
diffusion coefficient of water. The observed 
value of the broadening can be explained for a 
value of T,=4x10™* sec. Besides, the theo- 
retical quasi-elastic scattering in case (ii) has 
certain interesting features which are in general 
agreement with experiment. In part II of this 
paper, inelastic scattering (hindered transla- 
tions only) of cold neutrons has been calculated 
using two different models of water: (a) a gas 
model and (b) a Debye model; and the results 
have been compared with experiment. 

The general shape of both the quasi-elastic 
and inelastic scattering of cold neutrons and 
the magnitude of the diffusive broadening seem 
to support a quasi-crystalline model of water. 


*On leave of absence from the Atomic Energy Estab- 
lishment, Trombay, Bombay, India. 


ATOMIC MOTIONS IN WATER BY SCATTERING 
OF COLD NEUTRONS. D. J. Hughes,” H. Pal- 
evsky, W. Kley, and E. Tunkelo, Brookhaven 
National Laboratory, Upton, New York (Re- 
ceived March 29, 1960). 


The inelastic scattering of cold neutrons, of 
energy about 4x10~° ev, has been used to study 
the atomic motions in water, mainly in the 
liquid phase. As a result of the incoherent 
nature of scattering by protons, the interpreta- 
tion of the energy changes in terms of atomic 
motions is particularly simple. The water sam- 
ples used were extremely thin in order to avoid 
multiple scattering effects. Instead of the smooth 
distribution in energy of the scattered neutrons 
expected for a classical liquid, the experimental 
results exhibit a number of distinct energy 
changes. The observed transition energies are 
(in units of 10-° ev) 61, 21, 8, 5, and 0.7; the 
first three of these agree with Raman spectro- 
scopy results. The highest energy transition 
increases in intensity up to the boiling point, 
but does not shift in energy. For water vapor 
at 20 atmospheres, a smooth energy distribu- 
tion is obtained, which agrees well with the 
theory of Krieger and Nelkin. The sharp energy 
levels observed for liquid water indicate that the 
molecules do not act as free gas atoms of mass 
18; in the vapor the effective mass is about 4. 
An elastic peak is observed that exhibits no 
spread in energy of the type expected from the 
theory of classical diffusive motions. Instead, 
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two small but sharp peaks are found to corre- 
spond to the gain and loss of 0.7x10* ev energy; 
the nature of the transition corresponding to 
this energy change in unknown. The absence of 
diffusive broadening shows that the water mole- 
cules remain in one location for a relatively 
long time, about 10~'* sec, before undergoing 
diffusive “jumps.” In general, the present ex- 
periments show that the atomic motions in water 
are similar to those in a solid rather than a gas. 


‘Deceased. 


ELECTRICAL AND THERMAL CURRENTS IN 
A SLIGHTLY IONIZED GAS. Mahendra Singh 
Sodha, Physics Division, Armour Research 
Foundation, Chicago, Illinois (Received March 
31, 1960). 


In this communication, the author has solved 
the Boltzmann equation for electrons in a slightly 
ionized gas under the influence of an electric 
field (ac+dc), a magnetic field, and tempera- 
ture gradient. Expressions for the electrical 
and thermal currents have been obtained in 
terms of integrals having the collision frequency 
and/,, the isotropic part of the distribution 
function, in the integrand. A differential equa- 
tion for f, has been set up and analytical ex- 
pressions obtained under simplifying assump- 
tions. The application of the analysis to trans- 
port properties and electromagnetic wave 
propagation has also been indicated. 


ELECTROMAGNETIC SIGNALS FROM NUCLEAR 
EXPLOSIONS IN OUTER SPACE. Montgomery H. 
Johnson, Aeronutronic Division, Ford Motor 
Company, Newport Beach, California, and 
Bernard A. Lippmann, Lawrence Radiation Lab- 
oratory, University of California, Livermore, 
California (Received March 18, 1960). 


The thermal x rays produced by a nuclear burst 
in outer space cause polarization currents in the 
medium which, if distributed anisotropically, will 
emit electromagnetic radiation. Roughly, a burst 
of thermal x rays, equivalent in energy to 1 ton 
of high explosive, produces a detectable 10-Mc/ 
sec signal at a range of 1 km. Since only the 
ratio of x-ray energy to range enters into the 
strength of the radiated signal, other ranges fol- 
low by adjusting the x-ray energy proportionately. 
This works up to ~3 x10° km; beyond this range, 
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dispersive effects begin to reduce the signal re- 
ceived. 

The power in the electromagnetic signal varies 
as the square of the electron density, so this 
effect may provide a sensitive measure of the 
density of electrons in outer space. 





COMMENT ON A PAPER OF MORI ON TIME- 
CORRELATION EXPRESSIONS FOR TRANS- 
PORT PROPERTIES. Melville S. Green, 
National Bureau of Standards, Washington, D. C. 
(Received March 21, 1960). 





An autocorrelation expression given by Mori 
for the thermal conductivity of a fluid is shown 
to be only apparently different from an expres- 
sion previously derived by the author. 


CLUSTER INTEGRALS AND THE GROUND STAT: 
OF BOSONS WITH REPULSIVE INTERACTIONS. 
H. A. Gersch and V. H. Smith, School of Physics, 
Georgia Institute of Technology, Atlanta, Georgia 
(Received April 21, 1959; revised manuscript 
received April 19, 1960). 


The properties of a Bose system of particles 
with repulsive interactions have previously been 
treated using perturbation theory in the formalism 
of second quantization. Others have also consid- 
ered this problem by dealing with the wave func- 
tion in configuration space, using the theory of 
cluster expansions. In these latter papers, vari- 
ation with respect to a parameter in a trial func- 
tion for the ground state has been shown to yield 
a ground-state energy close to the exact asymp- 
totic expressions obtained from perturbation 
theory. The connection between the two methods 
is not immediately obvious from these cluster 
expansion treatments. It is shown here that, as 
one might expect, the cluster integral method can 
be handled so that it is completely equivalent to 
the pair approximation in perturbation theory. 


LATTICE DYNAMICS OF ALKALI HALIDE 
CRYSTALS. A. D. B. Woods, W. Cochran, and 
B. N. Brockhouse, Physics Division, Atomic 
Energy of Canada Limited, Chalk River, Ontar- 
io, Canada (Received March 11, 1960). 


The paper comprises theoretical and experi- 
mental studies of the lattice dynamics of alkali 
halides. A theory of the lattice dynamics of 
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jonic crystals is given, based on replacement 

of a polarizable ion by a model in which a rigid 
shell of electrons (taken to have zero mass) can 
move with respect to the massive ionic core. 

The dipolar approximation then makes the model 
exactly equivalent to a Born-von Kdrmdn crystal 
in which there are two “atoms” of differing 
charge at each lattice point, one of the “atoms” 
having zero mass. The model has been special- 
ized to the case of an alkali halide in which only 
one atom is polarizable, and computations of 
dispersion curves have been carried out for 
sodium iodide. We have determined the dis- 
persion v(q) relation of the lattice vibrations in 
the symmetric [001], [110], and [111] directions 
of sodium iodide at 110°K, by the methods of 
neutron spectrometry. The transverse acoustic, 
longitudinal acoustic, and transverse optic 
branches were determined completely with a 
probable error of about 3%. The dispersion 
relation for the longitudinal optic (LO) branch 
was determined for the [001] directions with less 
accuracy. Frequencies of some important pho- 
nons with their errors (units 10’* cps) are: 

TA (0,0, 1] 1.22+ 0.04, LA [0, 0, 1] 1.824 0.06, 
TA[$, 3, 3] 1.52+0.05, LA [, §, 5] 2.324 0.06, 

T0 [0, 0, 0] 3.6,+ 0.1, TO [0, 0, 1] 3.8,+0.1, 

10 [3, 3, 3] 3.5,+0.1. The neutron groups cor- 
responding to phonons of the LO branch were 
anomalously energy broadened, especially for 
phonons of long wavelength, suggesting a re- 
markably short lifetime for the phonons of this 
branch. 


NUCLEAR MAGNETIC RESONANCE IN COPPER 
ALLOYS— ELECTRON DISTRIBUTION AROUND 
SOLUTE ATOMS. T. J. Rowland, Union Carbide 
Metals Company, Niagara Falls, New York 
(Received March 18, 1960). 


The effects of the addition to copper of a wide 
variety of B subgroup elements on the nuclear 
magnetic resonance absorption of copper are 
described. The resonance amplitude, which 
udergoes a sharp reduction upon alloying, is 
of special interest; its dependence upon solute 
valence and size argues decisively in favor of 
conduction electron charge redistribution (va- 
lence effects) as the dominant source of the 


lectric field gradients surrounding these solutes. 


furthermore, these gradients are shown to 
decrease only about as 1/r* rather than expo- 
lentially as had been supposed. Using the pro- 





portional change in the lattice parameter of the 
solid solution as a measure of the local strains 
surrounding a solute atom, only slight correla- 
tions between local strains and resonance ampli- 
tude were found. It is concluded that the origin 
of electric gradients around multivalent solutes 
in copper is almost purely an effect of conduction 
electron distribution and that this distribution is 
not of the exponentially screened Coulomb charge 
type. The spatially oscillating charge distribution 
derived and recognized by Friedel and recently 
elaborated by Kohn and Vosko and Friedel and 
co-workers satisfactorily explains the observa- 
tions. 


THEORY OF NUCLEAR RESONANCE INTENSITY 
IN DILUTE ALLOYS. W. Kohn* and S. H. Vosko, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received March 18, 1960). 


Experiments of Bloembergen and Rowland have 
shown that the intensity of the nuclear resonance 
signa! in metallic Cu decreases rapidly when 
small quantities of other elements are alloyed 
with it. These results require that each solute 
atom produces significant electric field gradients 
in its vicinity, sometimes affecting as many as 
85 neighboring Cu nuclei. In this paper we show 
that field gradients of approximately the required 
magnitude arise from the redistribution of the 
conduction electron charge density near the solute 
atoms. A crucial feature of our theory is that 
at large distances y from a solute atom the 
electron density behaves as cos(2k°r+ @/r°, 
where k° is the Fermi wave number and ¢ is a 
phase. Our agreement with experiment is a 
confirmation of this behavior. Such an oscillatory 
behavior is a consequence of a discontinuous 
drop at the Fermi surface of n(k), the occupation 
probability of the conduction band function with 
wave vector k. 


* 
Present address: Department of Physics, University 
of California, La Jolla, California. 


PHONON SCATTERING IN KCI1-KBr SOLID 
SOLUTIONS AT LOW TEMPERATURES. 
Wendell S. Williams, * Department of Physics, 
Cornell University, Ithaca, New York (Re- 
ceived March 31, 1960). 


The scattering of phonons by point defects at 
liquid helium temperatures has been studied. 
Single crystals of KC1-KBr solid solutions were 
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employed since it is known that the Cl and Br™ 
ions are randomly distributed in the anion lattice 
sites. A 1 mole percent addition of KBr reduced 
the conductivity of KCl by a factor of three. For 
all the crystals measured, including pure KCl 
and pure KBr, the thermal conductivity showed 
a maximum at the same temperature (5°K). 
Hence, point defects reduce the conductivity on 
both the high- and low-temperature sides of the 
maximum. For the mixed crystals the thermal 
conductivity was found to be related to the ab- 
solute temperature by the empirical formula 

K =K,e -0.099T for the temperature interval 6°K 
to 16°K. 


"Now at Research Laboratory, National Carbon Com- 
pany, Division of Union Carbide Corporation, Parma, 
Ohio. 


DIFFUSION OF Li IN Si AT HIGH TEMPERA- 
TURES AND THE ISOTOPE EFFECT. E. M. 
Pell, General Electric Research Laboratory, 
Schenectady, New York (Received March 18, 
1960). 


The diffusion rate of Li in Si at high tempera- 
tures has been re-investigated using an out- 
diffusion technique. The resulting D for Li’ is 
(2.21+ 0.07) x10~* cm?/sec at (800+ 5)°C and 
2.4107 cm?/sec at (1350+5)°C. If these re- 
sults are combined with ion drift results, the 
diffusion constant can be described by 


D =(2.5+ 0.2) x107%exp| -(0.655+ 0.01)e/kT']. 


The isotopic effect upon the diffusion has been 
investigated using Li® and Li’. At 800°C, the 
value for Dy j0/Dy jn is 1.075+ 0.02, in accord- 
ance with the expected inverse dependence on 
the square root of the mass. The ionic charge 
of the Li and the atomic mechanism for Li dif- 
fusion are discussed in the light of these and 
other results. 


EXCITATION AND HIGH-TEMPERATURE AB- 
SORPTION OF KCI:Tl. David A. Patterson, 

U. S. Naval Research Laboratory, Washington, 
D. C. (Received March 22, 1960). 


Optical measurements have been made on crys- 
tals of KCl with a wide range of thallium concen- 
trations. Absorption measurements have been 
made up to 560°C and excitation spectra for 
luminescence have been measured from liquid 
nitrogen temperature to 100°C. These measure- 
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ments, in conjunction with earlier work, lead to 
the conclusion that there are at least seven bands 
in KC1:Tl: three in the “A” band region and two 
each in the “B” and “C” band regions. The large 
shift from “C” band to “B” band at high tempera- 
tures which has been previously reported was 
not found here. It is noted that the addition of 
small amounts of Sr to KCl shifts the fundamenta! 
absorption edge to short wavelengths. 


BROADENING OF SPIN- PHONON RESONANCE 
LINES BY EXCHANGE AND MAGNETIC DIPOLE 
INTERACTIONS. Rodney Loudon, Department 
of Physics, University of California, Berkeley, 
California (Received March 28, 1960). 


The second and fourth moments of the three 
principal ultrasonic free-spin absorption lines 
are calculated using a phenomenological form 
for the spin-phonon interaction. Both exchange 
and dipole interactions are taken into account, 
and it is found that exchange causes increased 
linewidth in all three cases. For the line at the 
Larmor frequency, the moments are compared 
with those for the corresponding photon absorp- 
tion line, for which exchange narrowing occurs. 


AUGER ELECTRON EJECTION FROM GER- 
MANIUM AND SILICON BY NOBLE GAS IONS. 
Homer D. Hagstrum, Bell Telephone Labora- 
tories, Murray Hill, New Jersey (Received 
March 21, 1960). 


Experimental results concerning electron 
ejection from annealed, atomically clean sur- 
faces of germanium and silicon by the singly 
charged ions of the noble gases are reported. 
The (111) and (100) faces of silicon and the (111) 
face of germanium have been studied. Total 
yield and kinetic energy distribution of ejected 
electrons were measured and ion energies varied 
in the range 10 to 1000 ev. A new method of 
operation of the apparatus and of obtaining the 
kinetic energy distributions from the recorded 
retarding potential data has been employed. 
Documentation of the state of the target surfaces 
is given including photomicrographs and electron 
micrographs of the silicon surfaces. Since these 
experimental results are subsequently to be 
interpreted theoretically, identification of the 
results with the theoretical ideas only is given 
here. 
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PARAMAGNETIC RESONANCE OF V** IN TiO,. 
H. J. Gerritsen and H. R. Lewis, RCA Labora- 

tories, Princeton, New Jersey (Received March 
28, 1960). 


The paramagnetic spectrum of vanadium in TiO, 
has been observed at 10.14 and 22.68 kMc/sec. 
An analysis of the data indicates that the spec- 
trum is due to single d electrons of tetravalent 
vanadium ions located at titania sites in the 
lattice. The calculated values of the components 
of the g tensor and the hyperfine interaction con- 
stant are: g, =1.915, gy =1.912°, g, =1.956°, 

A, =0.0031° cm™, Ay =0.0043 cm™, and Az 
=0.0142 cm~. The axes of the magnetic co- 
ordinate systems of the two nonequivalent ions 
per unit cell are in the [1, 1, 0] and [ 1, 1, 0] 
directions. 


QUENCHING OF LATTICE VACANCIES IN PURE 
SILVER. M. Doyama and J. S. Koehler, Univer- 
sity of Illinois, Urbana, Illinois (Received March 
30, 1960). 


Quenching data on 99.999% pure silver give 
the energy required to form a lattice vacancy to 
be 1.10+ 0.04 ev. In addition, the data of Simmons 
and Balluffi who obtained the vacancy concentra- 
tion in equilibrium near the melting temperature 
together with present data give the resistivity 
increase per atomic percent vacancies to be 1.3 
+0.7 micro-ohm cm. 


ELECTRON SPIN-LATTICE RELAXATION IN 
DILUTE POTASSIUM CHROMICYANIDE AT 
HELIUM TEMPERATURES. J. G. Castle, Jr., 

P. F. Chester, and P. E. Wagner,* Westing- 
house Research Laboratories, Pittsburgh, Penn- 
sylvania (Received March 18, 1960). 


Measurements have been made of the electron 
spin-lattice relaxation of the -1/2, +1/2 line of 
Cr*** in K,Co(CN), at 9 kMc/sec as a function 
of temperature, chromium concentration, and 
the proximity of the -3/2, -1/2 line. The ex- 
perimental procedure, involving inversion of the 
line, is capable of distinguishing a “bottleneck” 
relaxation time from a true spin-phonon relaxa- 
tion time, T,. At Cr*** concentrations up to 
0.5%, the relaxation data are fitted well by single 
exponential functions of time. Between 1.3°K and 
4.8°K, T, varies approximately as T~*? indicating 
that the single-phonon process is dominant. No 





phonon-bath bottleneck is observed, in agreement 
with calculations based on the measured param- 
eters. A “proximity effect” is observed in which 
the relaxation rate of the -1/2, +1/2 line is en- 
hanced when the -3/2, -1/2 line is within 20 
linewidths. At 1% Cr***, the relaxation behavior 
is markedly different: the recovery is consider- 
ably faster and can no longer be described by a 
single time constant. This change and the prox- 
imity effect are interpreted qualitatively in 
terms of spin cross relaxation. The measured 
linewidth increases with concentration from 
0.03% to 2% Cr***, even though the line is ob- 


served to be inhomogeneous at and below 0.5% 
Crttt, 


*Present address: Department of Electrical En- 
gineering, The Johns Hopkins University, Baltimore, 
Maryland: 


ELECTRON WAVE FUNCTIONS IN METALLIC 
POTASSIUM. Joseph Callaway, Department of 
Physics, University of Miami, Coral Gables, 
Florida (Received March 18, 1960). 


Wave functions to order k? are presented for 
electrons in metallic potassium. The calculation 
is an application of the cellular method. The 
potential was derived from a self-consistent 
field and contains exchange effects. 


TRIGONAL MAGNETOCRYSTALLINE ANISOT - 
ROPY IN HEXAGONAL OXIDES. L. R. Bick- 
ford, Jr., International Business Machines, 
Yorktown Heights, New York (Received March 
29, 1960). 


Torque measurements of the magnetocrystalline 
anisotropy between 77°K and 300°K are reported 
for single crystals of two different ferrimagnetic 
oxides having structures related to that of the 
mineral magnetoplumbite. The compounds, 
with chemical compositions Co,Ba,Fe,,0,, and 
Co,Ba,Fe,,0,,, are known as Co,Y and Co,Z, 
respectively. Both compounds, after suitable 
magnetic field cooling treatment, display trigon- 
al anisotropy in the basal plane at 77°K. The 
field cooling is shown to have the effect of plac- 
ing the magnetization into one of two energetically 
equivalent orientations, each of which leads toa 
trigonal term of different sign. Examination of 
the crystal structures shows that in the case of 
Co,Y an additional term K,’sin*6cos@cos3y 
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(where @ and ¢ are the polar coordinates of the 
magnetization) should be added to the generally 
accepted phenomenological hexagonal anisotropy 
expression. The anisotropy constant K,’ is 
evaluated as 610° ergs/cm® at 117°K. In the 
case of Co,Z this trigonal anisotropy term is 
not consistent with the crystal symmetry, al- 
though it is appropriate for structural sub- 
groups of the unit cell. Its trigonal anisotropy 
is explained in terms of a “puckered” magneti- 
zation pattern whereby the sign of the c-axis 
component of magnetization is different for 
adjacent subgroups. The existence of this puck- 
ered pattern implies that the exchange coupling 
across the boundaries between subgroups is 
relatively weak. A new rotational hysteresis 
effect in Co,Y is described and explained pheno- 
menologically. An atomic theory assigning the 
origin of the trigonal anisotropy of both com- 
pounds to the cobalt ions is presented. 


GALVANOMAGNETIC AND THERMOMAGNETIC 
POTENTIALS IN ZINC AT LIQUID HELIUM 
TEMPERATURES. C. J. Bergeron, C. G. 
Grenier, and J. M. Reynolds, Department of 
Physics, Louisiana State University, Baton 
Rouge, Louisiana (Received January 21, 1960; 
revised manuscript received April 21, 1960). 


The Hall effect, magnetoresistance, thermo- 
electric voltage, and transverse Ettinghausen- 
Nernst potential have been measured in a single 
crystal of zinc in the liquid helium temperature 
range. Oscillations as a function of magnetic 
field strength were observed in all of these 
potentials. The measured period in 1/H for the 
transverse oscillations was 6.2 107° gauss 
+0.5%, with the magnetic field parallel to the 
hexagonal axis. Both transverse effects pos- 
sessed strong second harmonic oscillations. 
The oscillations in the longitudinal effects both 
exhibited a phase inversion in the neighborhood 
of 4.2 kilogauss, the same field region for which 
the gross Hall field changed sign. In this same 
field region the period of the oscillations for the 
longitudinal effects was that of the second har- 
monic. 

Empirical correlations between the reversal 
of sign of the Hall effect and (1) the phase re- 
versal of the magnetoresistance oscillations, 

(2) the strong second harmonic content of these 
oscillations in the region of the phase reversal, 
and (3) the quadratic shape of the envelope of 
the magnetoresistance oscillations in this region 
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can be achieved by assuming the oscillatory 
component of o,, to be much more significant 
than that of o,,, where o;, is the conductivity 
tensor. 

A further correlation between the oscillatory 
thermal effects can be achieved by assuming 
that the difference of the absolute thermoelectric 
power and the temperature derivative of the 
chemical potential is negligible. 







EFFECT OF UNIFORM COMPRESSION ON THE 
ETTINGHAUSEN-NERNST EFFECT IN ZINC AT 
LOW TEMPERATURES. K. S. Balain, C. G. 
Grenier, and J. M. Reynolds, Department of 
Physics, Louisiana State University, Baton 
Rouge, Louisiana (Received January 29, 1960). 





The influence of hydrostatic pressure on the 
oscillatory Ettinghausen-Nernst effect in a zinc 
single crystal has been studied in magnetic 
fields ranging from 2 to 11 kilogauss at 4.2°K. 
With the field along the hexagonal axis, the long- 
period oscillations decrease in period by 5% as 
the pressure is increased from 300 to 2300 psi. 
There is an average increase in amplitude of 
about 10% for each 250 psi increase of pressure. 
There is no evidence of phase change with pres- 
sure. 


SPECTRAL DISTRIBUTION OF THE PHOTO- 
MAGNETOELECTRIC EFFECT IN Ge: EXPER- 
IMENT. Frank A. Brand, Andrew N. Baker,* 
and Herbert Mette, U. S. Army Signal Research 
and Development Laboratory, Fort Monmouth, 
New Jersey (Received March 31, 1960). 


The photomagnetoelectric effect has been 
studied in germanium as a function of the wave- 
length of incident radiation in the region from 
0.5 to 2.0 microns. The dependence of both 
photoconductivity and photomagnetoelectric 
(PME) response has been measured in various 
samples, using front and back surface recom- 
bination velocities and bulk recombination as 
parameters. It has been found that under cer- 
tain conditions of bulk and surface recombination 
a reversal in sign for the PME response occurs 
over the frequency range studied. Sign rever- 
sals obtain at wavelengths in the range from 1.55 
to 1.85 microns, corresponding to partial optical 
transparency. The exact frequencies at which 
reversal occurs depend on the surface and bulk 
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recombination rates, the condition being that 

the Dember field be zero corresponding to equal 
carrier concentrations on the front and back 
surfaces. Furthermore, the present observa- 
tions are shown to be in good qualitative agree- 
ment with the theoretical work reported by 
Gartner. Experimental procedures are described, 
and it is shown how this effect can be used to 
advantage in the study of surface recombination 
velocities in various environments. 


*Now at Bell Telephone Laboratories, Murray Hill, 
New Jersey. 


THEORY OF SOLID He*. N. Bernardes, Insti- 
tute for Atomic Research and Department of 
Physics, Iowa State University, Ames, Iowa, 

and H. Primakoff, Department of Physics, Wash- 
ington University, St. Louis, Missouri (Received 
March 30, 1960). 


A theoretical analysis is given of the properties 
of solid He*® on the basis of: (1) a gas-phase 
Lennard-Jones “12-6” potential modified at small 
interatomic distances; (2) a Heitler- London type 
variational-trial wave function for all the atoms 
in the solid constructed from a properly anti- 
symmetrized product of individual atom orbitals 
localized on the various lattice points; (3) a 
Dirac vector model to describe the symmetry 
energy with an exchange integral deduced from 
(1) and (2); (4) a spin wave approximation at 
“low” temperatures and a Kramers-Opechowski 
approximation at “high” temperatures for calcu- 
lation of the free energy of the nuclear spins; 
and (5) a Debye phonon model for the description 
of the vibrationally excited states of the solid. 

On this basis, calculated values at low pressures 
and temperatures (p~=30 atm; T$1°K) are pre- 
sented for: (a) the cohesive energy per atom: 
=2.5°K; (b) the root mean square deviation of an 
atom from its lattice site: =0.36 x nearest neigh- 
bor distance; (c) the nuclear magnetic suscepti- 
bility which corresponds to an antiferromagnetic 
behavior with a “paramagnetic” Curie tempera- 
ture T, =0.1°K; (d) the variation (decrease) of 

Tc with increasing pressure corresponding to a 
possible nuclear antiferromagnetic to nuclear 
ferromagnetic transition for p~=150 atm; (e) the 
Specific heat which exhibits an anomaly at 
T=0.1°K associated with the alignment of the nu- 
clear spins; (f) the thermal expansion coefficient 
which becomes negative below about 0.6°K; (g) 

the melting curve which is characterized by a 


minimum at T=0.37°K and a maximum at 
T =0.08°K. Comparison of the theory is made 
with available experimental data. 


PRESSURE BROADENING OF THE VIOLET 
TRIPLET OF MANGANESE IN THE PRESENCE 
OF ARGON AND HELIUM. Shang-Yi Ch’en and 
Robert B. Bennett, * University of Oregon, Eu- 
gene, Oregon (Received June 29, 1959; revised 
manuscript received March 21, 1960). 


The 4030 triplet of Mn (3d°4s?, °S - 3d54s4p, °P°) 
was studied in absorption under various pres- 
sures of argon and helium ranging in relative 
density (r.d.) from 1 to 30. Helium produced no 
shift while argon produced a linear red shift of 
0.12 cm~/r.d. In the presence of both helium 
and argon the half-widths were found to increase 
directly with increasing relative density of the 
perturbing gas. The j=3/2, 5/2, and 7/2 com- 
ponents are broadened in helium by 0.46, 0.50, 
and 0.47 cm™/r.d., respectively, while the 
same components are broadened in argon by 
0.33, 0.36, and 0.30 cm/r.d. The result for 
Mn/Ar indicates that the interatomic energies 
are proportional to the inverse sixth power of 
the interatomic distance and that the difference 
in the interaction constants for the ground state 
®S and the excited states °P of Mn in argon is 
3x10 erg cm*. 


*Now at Whitman College, Walla Walla, Washington. 


HYPERFINE STRUCTURE OF THE MICRO- 
WAVE SPECTRA OF THE NO MOLECULE AND 
THE NUCLEAR QUADRUPOLE MOMENT OF 
NITROGEN. Chun C. Lin, Department of Phys- 
ics, University of Oklahoma, Norman, Oklahoma 
(Received March 23, 1960). 


The frequencies of the magnetic resonance spec- 
trum of the NO molecule have been recalculated 
by using the new values of the spin-orbit coup- 
ling constants and by taking the effect of 7 un- 
coupling into consideration. The agreement be- 
tween the theoretical and experimental results 
is improved over the previous calculation. By 
combining the magnetic hyperfine and nuclear 
quadrupole coupling constants, the ratio of the 
quadrupole moment to the magnetic moment of 
the nitrogen nucleus is obtained. The nuclear 
quadrupole moment of nitrogen is found to be 
(0.016 + 0.007) x10-** cm*. The uncertainty of 


593 








VoLuME 4, NuMBER 11 





PHYSICAL REVIEW LETTERS 





JuNE 1, 1960 





this value is chiefly due to that of the coupling 
constants rather than to the nature of the method 
itself. 


VIBRATIONAL STATES OF THE HYDROGEN 
MOLECULAR ION. Stanley Cohen, John R. 
Hiskes, and Robert J. Riddell, Jr., Lawrence 
Radiation Laboratory, University of California, 
Berkeley, California (Received March 31, 1960). 


The eigenvalues and eigenfunctions of the vi- 
brational states belonging to the ground elec- 
tronic state of the hydrogen molecular ion have 
been calculated. The calculations have been 
done for the J=0, 2, 4, and 7 rotational states. 
Included is a discussion of the dependence of the 
eigenvalues as a function of the lowest-order 
dynamic corrections to the internuclear poten- 
tial. Also, a calculation has been done to deter- 
mine the number of bound states of the D,* sys- 
tem. 


NONLOCAL OPTICAL MODEL FOR NUCLEON- 
NUCLEAR INTERACTIONS. P. J. Wyatt, Aero- 
nutronic, Ford Motor Company, Newport Beach, 
California, J. G. Wills, Los Alamos Scientific 
Laboratory, Los Alamos, New Mexico, and 

A. E. S. Green, Physics Section, Convair -San 
Diego, San Diego, California (Received Decem- 
ber 8, 1959). 


An attempt is made to achieve a unified poten- 
tial description of the gross structure of the 
nucleon-nuclear interactions in bound states and 
in states of scattering. A model is employed 
with a nonlocal complex diffuse potential with 
spin-orbit coupling and surface absorption. This 
represents a relatively simple nonlocal generali- 
zation of the usual static models which might 
reasonably be expected to describe the nucleon- 
nuclear interaction in the low-energy range (say 
from -25 Mev to 25 Mev). Choosing the range 
of the nonlocal forces as suggested by considera- 
tions of the properties of infinite nuclear matter, 
the real parameters are fixed largely on the 
basis of neutron and proton separation energies. 
Two absorption parameters are then adjusted to 
total reaction and differential elastic cross- 
section data for neutrons. It is found that the 
successes of local optical models with energy- 
dependent parameters are largely preserved 
despite the rather large range of nonlocality. 
Contrary to expectations, it is found that nonlo- 
cality tends to accentuate rather than wash out 
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diffraction patterns. Although a diverse variety 
of experimental phenomena are treated, a range 
of parameter choices remains. Because of 
theoretical uncertainties as to the size of the 
“rearrangement energy,” an effort is made to 
establish limits as to its magnitude on pheno- 
menological grounds. The influence of several 
choices upon the physical phenomena used in 
adjusting the parameters of this model is shown. 
It would appear that this study does allow for a 
rearrangement energy but that it is rather small 
(< 6 Mev) and comparable to the fluctuations of 
the potential from element to element. 


SINGLE-PARTICLE STATES IN DEFORMED 
NONLOCAL DIFFUSE BOUNDARY POTENTIALS. 
R. H. Lemmer,* Physics Department, Florida 
State University, Tallahassee, Florida, and 

A. E. S. Green, Physics Section, Convair-San 
Diego, San Diego, California (Received February 
24, 1960). 


Using the spherical wave functions generated 
in a previous investigation by Wyatt, Wills, and 
Green, the influence of spheroidal deformation 
is examined with the aid of perturbation theory. 
The combined calculations yield the energies of 
single-particle states for a diffuse-boundary, 
nonlocal deformed potential. Specific calcula- 
tions are performed for light nuclei around A =25 
and in the rare earth region between A = 150 and 
A=180. An analysis of nuclear ground-state 
spins and magnetic moments is presented in 
terms of the computed level schemes and wave 
functions. The results confirm the general as- 
pects of the Nilsson, Mottelson results as ob- 
tained with adjusted harmonic oscillator poten- 
tials although some differences arise in detail. 
In particular, the calculated coefficients usually 
show less mixing of different angular momentum 
states in our case. The fact that the unperturbed 
potentials used in this calculation were obtained 
in the study of Wyatt, Wills, and Green from 
completely independent theoretical and experi- 
mental considerations is satisfying and further 
tends to confirm that the phenomenological model 
has a strong basis in reality. A discussion of 
the relationship of the phenomenological model 
to the self-consistent nuclear model of Brueckner 
is given. 


* 
Now at Physics Department, Massachusetts Insti- 
tute of Technology, Cambridge, Massachusetts. 
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EXCITATION CURVES AND ANGULAR DISTRI- 
BUTIONS FOR N"*(d,n)O'5. Theo Retz-Schmidt 
and Jesse L. Weil, The Rice Institute, Houston, 
Texas (Received March 24, 1960). 


Excitation curves for the highest energy neu- 
tron group in the reaction N“*(d,”)O"* have been 
measured at 6),,,=0°, 30°, 90°, and 164° for 
deuteron bombarding energies between 0.66 and 
5.62 Mev. A pulse shape discrimination detector 
was used to eliminate the pulses due to y rays 
from the neutron spectra. There is considerable 
resonance structure in the excitation curves, 
with the anomalies appearing at different en- 
ergies for the different angles. The angular 
distribution of this neutron group has also been 
measured at bombarding energies of 0.91, 1.17, 
1.51, 1.88, 2.58, 3.13, 3.56, 4.36, 4.80, and 
5.27 Mev. The shape of the angular distribution 
changes rapidly with energy at low bombarding 
energy, but above 3.5 Mev the shape becomes 
more stable. The maximum cross section at 
any angle was 5.5 millibarns per steradian. 


EFFECTS OF CHEMICAL BINDING ON NU- 
CLEAR RECOIL. M. S. Nelkin and D. E. Parks, 
John Jay Hopkins Laboratory for Pure and 
Applied Science, General Atomic, San Diego, 
California (Received March 28, 1960). 


The recoil of a chemically bound nucleus is 
considered for slow neutron scattering and for 
the resonant absorption of neutrons or gamma 
rays. The Doppler-broadened resonance line 
shape is derived in terms of the time-dependent 
self-correlation function describing the motion 
of a nucleus due to the interatomic forces. This 
explicitly relates the resonance line shape to 
the differential scattering cross section for slow 
neutrons in the Fermi pseudopotential approxi - 
mation. Within this formulation an expansion 
for large nuclear recoil is naturally suggested. 
For the case of a crystal, this expansion can 
be directly related to the expansion associated 
with the central limit theorem of probability 
theory and can therefore be proved to be asymp- 
totic in nature. The expansion parameter is 
(K,,/R)”, where Kay is the average kinetic 
energy of a nucleus and R is the recoil energy 
for a free nucleus at rest. The leading term of 
the expansion is the weak-binding limit origin- 
ally obtained by Lamb. In this limit the Doppler - 
broadened line shape is the same as would ob- 
tain for an ideal monatomic gas of the same 





mass with an effective temperature T’ =(2/3 )\Kay- 
For noncrystalline systems, a similar expan- 
sion with the same leading term can be obtained 
by a rearrangement of the terms in an expansion 
used by Wick to study the slow neutron total 
cross section. The relation of the present ex- 
pansion to Wick’s expansion is discussed. 


HYPERFINE STRUCTURE AND NUCLEAR 
MOMENTS OF 17-hr BROMINE-76. Edgar Lip- 
worth, Thomas M. Green, Hugh L. Garvin, and 
William A. Nierenberg, Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received March 25, 1960). 


The nuclear spin, the nuclear magnetic -dipole 
interaction constant a, and the nuclear electric- 
quadrupole interaction constant b have been de- 
termined for 17-hr bromine-76 by an atomic- 
beam experiment. The results are: J=1, |a| 
=345.422+ 0.014 Mc/sec, |b| =314.329+ 0.022 
Mc/sec, 6/a=0.9100+ 0.0001. The nuclear mag- 
netic-dipole and electric-quadrupole moments 


are calculated to be, respectively, » =+0.5479 
+ 0.0001 nuclear magneton and @ =¥0.27+0.01 


barn. The sign of u, though not determined, 
is probably negative. 

The hyperfine structure separations are 
Av(5/2, 3/2) =1256.47+ 0.05 Mc/sec and 
Av(1/2, 3/2) =189.11+0.05 Mc/sec. The hyper - 
fine structure is of particular interest because 
the F =1/2 and 3/2 levels are inverted and not 
in normal order. This inversion is the first 
case of its kind established in an atomic-beam 
experiment. 


ANGULAR DISTRIBUTIONS OF B’°(d, a)Be® RE- 
ACTIONS FROM 0.6 TO 1.5 Mev. Robert L. 
Becker,* University of Kentucky, Lexington, 
Kentucky (Received December 7, 1959). 


Angular distributions of the two most energetic 
alpha-particle groups resulting from deuteron 
bombardment of B*® were measured at seven deu- 
teron energies between 0.58 and 1.50 Mev. In 
addition to the usual pulse-height analysis of the 
detector output, pulse-decay analysis was also 
employed so that alpha particles could be dis- 
tinguished from protons giving the same pulse 
size. The shapes of the angular distributions 
were somewhat energy dependent. Averaged 
over the range of energies studied, the angular 
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distributions were not isotropic, the yield being 
slightly higher at back angles. 


“Now at Boston College, Chestnut Hill, Massachu- 
setts. 


CALCULATION OF a-TRANSITION PROBABIL- 
ITIES. Hans J. Mang, Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received February 15, 1960). 


The decay rates for the ground-state transi- 
tions of all polonium isotopes and the odd-even 
astatine isotopes are discussed on the basis of 
the nuclear shell model. Good agreement with 
experimental data is obtained. In particular the 
behavior of the reduced width as a function of 
the neutron number around the magic number 
N =126 is well reproduced. 


* 
Present address: University of Heidelberg, Heidel- 
berg, Germany. 


K -DEUTERON SCATTERING AND THE K’ - 
NUCLEON SCATTERING LENGTHS. T. B. Day, 
G. A. Snow, and J. Sucher, University of Mary- 
land, College Park, Maryland (Received March 
16, 1960). 


Cross sections for K -d reactions have been 
calculated in the low-momentum region for sev- 
eral possible values of the elementary K -nucleon 
scattering amplitudes. Multiple-scattering ef- 
fects have been included in an approximate way. 
A comparison of the results for the sum of the 
elastic plus breakup cross sections with the pre- 
liminary measurements available is presented. 


CHARGE -EXCHANGE SCATTERING OF NEGA- 
TIVE PIONS AT 150 Mev. W. J. Kernan,” The 
Enrico Fermi Institute for Nuclear Studies and 
The Physics Department, The University of 
Chicago, Chicago, Illinois (Received March 31, 
1960). 


The charge-exchange scattering of negative 
pions by hydrogen has been measured at a bom- 
barding energy of 150 Mev. The energy dis- 
tribution of gamma rays from the decay of the 
neutral pions was measured with a lead glass 
Cerenkov counter at laboratory angles of 45°, 
75°, 105°, and 135°. If the charge -exchange dif- 
ferential scattering cross section in the center- 
of-mass system is expanded as a series of 
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Legendre polynomials, the result is 
do/dQ =(1.00+ 0.03)[3.39+ 0.11 
- (1.54 + 0.29)P,(cosé@’) 
+(3.57+ 0.56)P,(cosé’)} mb-sr~. 


The total cross section for charge exchange, 
obtained by integration, is then o%,4(7 —7°) 
=42.6+1.9 mb. 


"Present address: Argonne National Laboratory, 
Lemont, Illinois. 


CHARGE -EXCHANGE SCATTERING OF NEGA- 
TIVE PIONS AT 61 Mev AND 95 Mev. C. M. 
York,* W. J. Kernan, and E. L. Garwin,! The 
Enrico Fermi Institute for Nuclear Studies 

and The Physics Department, The University of 
Chicago, Chicago, Illinois (Received March 31, 
1960). 


The charge-exchange scattering of negative 
pions by liquid hydrogen has been measured at 
61+ 1 Mev and 95+ 2 Mev bombarding energy. 
The measurements were made with a gamma- 
ray spectrometer which employs a lead glass 
Cerenkov counter. If the charge-exchange scat- 
tering cross section is expanded as a series of 
Legendre polynomials in the center-of-mass 
system of the collision, we find that at 61 Mev, 


do/dQ =(1.00+ 0.05)[0.613 + 0.030 

- (0.830 + 0.068)P, (cos@’) 

+ (0.183 + 0.150)P,(cos@’)] mb-sr™, 
and at 95 Mev, 
do /dQ = (1.00 + 0.03)[1.05 + 0.05 

- (1.15 + 0.12)P,(cosé@’) 

+ (0.33 + 0.25)P,(cos@’)] mb-sr™. 





The total cross section for charge exchange, 
obtained by integration, is: o4,4(7 —7°)=7.7+0.6 
mb at 61 Mev and o,4(7" ~ 7°) =13.2+ 0.8 mb at 
95 Mev. 

A table summarizing the measurements per- 
formed by this group at 61 Mev, 95 Mev, 128 
Mev, and 150 Mev is given. 


* 

Present address: CERN, Geneva, Switzerland. 

Present address: Argonne National Laboratory, 
Lemont, Illinois. 

tPresent address: Department of Physics, Univer- 
sity of Illinois, Urbana, Illinois. 
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ARRIVAL DIRECTION OF COSMIC -RAY AIR 
SHOWERS FROM THE EQUATORIAL SKY. E. V. 
Chitnis* and V. A. Sarabhai, Physical Research 
Laboratory, Ahmedabad, India, and G. Clark, 
Laboratory for Nuclear Science, Massachusetts 
Institute of Technology, Cambridge, Massachu- 
setts (Received March 22, 1960). 


The celestial arrival directions of over 100000 
showers with sizes greater than 10° particles 
have been determined by fast timing in observa - 
tions at an altitude of 2034 m. The observations 
covered a band of declinations from -30° to +50° 
with an angular resolution of 4°, and they ex- 
tended a survey begun in an earlier experiment 
that covered the northern sky. As in the earlier 
experiment no significant deviation from isotropy 
was found. The atmospheric attenuation of the 
shower intensity was determined from the zenith 
angle distribution, and also from a comparison 
of the absolute shower intensity at 2034 m and at 
sea level. Within an experimental uncertainty 
of about 5%, both methods yield an exponential 
attenuation length consistent with the value of 
107 g cm™~ previously found at sea level. The 
absolute intensity of showers with more than 10° 
particles at 2034 m was found to be (1.11+ 0.30) 
x10? cm™ sec™ sr“. 


‘Present address: Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cambridge, 
Massachusetts. 


PION- PION SCATTERING AND K*— 37 DECAY. 
N. N. Khuri,* Institute for Advanced Study, 
Princeton, New Jersey, and S. B. Treiman, 
Palmer Physical Laboratory, Princeton Univer- 
sity, Princeton, New Jersey (Received March 
16, 1960). 


The effects of final-state pion-pion interac- 
tions on the spectrum of K*~ 37 decay is studied 
by dispersion relation methods. In the approxi- 
mations adopted we are led to a set of linear in- 
tegral equations for the amplitudes of the K*—3z 
decay. The kernels in these equations depend on 
the pion- pion S-wave scattering amplitudes. An 
approximate solution for these equations is ob- 
tained by iteration and the departures from a 
purely statistical spectrum for the decay are re- 
lated to pion-pion S-wave scattering. The latter 
inturn is assumed to be well represented with a 
scattering length structure. The K*~ 37 spec- 
trum then is parametrized by two quantities, the 








T=0 and T=2 pion-pion S-wave scattering 









lengths, a, anda,. Such experimental results as 
presently exist indicate that a, - a, is positive 
and that roughly a, -a,~0.7, in units of the pion 
Compton wavelength. 


“On leave of absence from the American University 
of Beirut, Beirut, Lebanon. 


GAUGE INVARIANCE AND THE LORENTZ PON- 
DERMOTIVE FORCE. Lloyd Motz, Rutherfurd 
Observatory, Columbia University, New York, 
New York (Received March 3, 1958; revised 
manuscript received January 28, 1959). 


In this paper it is shown that if one introduces 
into the Weyl theory of gauge invariance the two 
additional conditions that gauge (and therefore 
length), except for an arbitrary phase factor, be 
integrable along the path of a particle, and that 
the change in dimensions of a particle be a mini- 
mum, one immediately obtains the Lorentz pon- 
dermotive force for a charged particle in an 
electromagnetic field and the Bohr-Sommerfeld 
quantum integral. 


FINAL-STATE INTERACTIONS AND |A1/| = 1/2 
SELECTION RULE. K. Chadan* and S. Oneda,t 
Department of Physics, University of Maryland, 
College Park, Maryland (Received March 28, 
1960). 


We consider the effect of various isotopic spin 
selection rules and final-state interactions be- 
tween two outgoing pions upon the K,4/K,3 
branching ratio. It is shown that in the leptonic 
decay modes of the K meson the transition am- 
plitude contains at most a small amount of | Al! 
=3/2. 


“On leave of absence from Ecole Normale Supérieure, 
Paris, France. 
TNow at University of Kanazawa, Kanazawa, Japan. 


UNSTABLE PARTICLES IN A GENERAL FIELD 
THEORY. J. Gunson and J. G. Taylor, Depart- 
ment of Applied Mathematics and Theoretical 
Physics, University of Cambridge, Cambridge, 
England (Received March 14, 1960). 


The problem of unstable particles in quantum 
field theory is treated as one of the interpreta- 
tion of complex singularities appearing in the 
analytic continuation of scattering amplitudes 
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into unphysical sheets of their Lorentz-invariant 
variables. Suitable continuations are shown to 
hold under certain restrictive assumptions in a 
general field theory, making use of unitarity and 
causality of the S matrix. The extra singulari- 
ties appearing in the continuation are fixed iso- 
lated poles, in accordance with a conjecture of 
Peierls. 


ELECTROMAGNETIC FORM FACTORS OF THE 
NUCLEON. F. J. Ernst, University of Wisconsin, 
Madison, Wisconsin and Columbia University, 
New York, New York, R. G. Sachs, University 

of Wisconsin, Madison, Wisconsin and Univer- 
sity of Paris, Paris, France, and K. C. Wali, 
University of Wisconsin, Madison, Wisconsin 
(Received March 18, 1960). 


The physical interpretation of the electromag- 
netic form factors is discussed with special re- 
ference to the gauge invariance of particular 
theories. A distinction is made between the 
condition that the one-nucleon matrix element 
satisfy the equation of continuity (“weak gauge 
invariance”) and the stronger condition imposed 
by the generalized Ward identity (“ strong gauge 
invariance”). The former is shown to be a con- 
sequence of covariance under the improper 
Lorentz transformations, and hence it has no 
new content concerning the functional behavior 
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of the form factors. The latter implies restric- 
tions on the current operator which may have an 
important effect on the results of calculations of 
form factors. 

In connection with the physical interpretation 
it is noted that the moments of the charge and 
current distribution are determined by F.,=F, 
-(q?/2M) F, and F = (1/2M) F,+F,. Specific- 
ally the second moment of the charge distribu- 
tion, - 6 F .,'(0), is found, in the case of the 
neutron, to be directly measured by the neutron- 
electron interaction without the intervening sub- 
traction of the Foldy term. 

These matters are investigated in detail by 
means of a specific model of the nucleon which 
is a covariant generalization of the fixed-source 
static model having the property that it gives 
results identical with the static model in the 
limit M-o. It is found that strong gauge invar- 
iance requires the addition of line currents 
which make significant contributions to the form 
factors in general and, in particular, to the 
proton charge radius even in the static approxi- 
mation. This suggests that as a consequence of 
strong gauge invariance, important contributions 
to the charge radius must arise in any theory 
from intermediate states of large mass. The 
model also provides a means of consistently 
calculating recoil corrections to the static model. 
They are found to be large. 











